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FREE VIBRATION OF A ROTATING 
BEAM-CONNECTED SPACE STATION* 


SUMMARY 


The free vibration of a rotating beam-connected space station is 
analyzed with a mathematical model of the space station which represents the 
general three-dimensional motion of tiie various components of the system. 

The space station is composed of two space modules connected by a flexible 
beam, and the system is caused to spin in the plane of its orbit in order to 
produce an artificial gravity environment within the space modules. 

The kinetic energy and potential energy of the space station are used to 
develop a Lagrangian function of the system. Hamilton's principle is used to 
determine a set of governing equations, and a set of boundary conditions repre- 
senting a clamped-clamped attachment of the beam to each space module is 
applied to the ends of the beam. Within the limits of small deflection theory, 
the motion of the space station is shown to be uncoupled into two separate 
types of motion, one in the plane of rotation and the other perpendicular to 
the plane of rotation. 

An exact solution is obtained for the beam deflection in the plane of 
rotation. The application of the nonhomogeneous boundary conditions leads 
to a set of simultaneous equations in the frequency p 2 , from which a charac- 
teristic determinant is developed. A procedure to solve for the zeros of the 
characteristic determinant is programmed for digital solution on the IBM 7094. 

Results of the analysis for a given space station design are presented 
in the form of tables showing the natural frequencies of free vibration of the 
space station for various spin rates. The effect of the spin rate is shown to 
be an "apparent" increase in the stiffness of the beam. Mode shapes showing 
the normalized bending deflection of the beam in its six lowest modes of 
vibration are presented. By analysis of five special cases of negative values 
of p 2 , the existence of two rigid body modes with nonzero values of p 2 is 
demonstrated; and it is shown that the configuration studied has no instabilities 
for motion in the plane of the orbit. 


* The information presented herein was included in a dissertation 
entitled, "Free Vibration of a Rotating Beam-Connected Space Station" sub- 
mitted in partial fulfillment of the requirements for the degree of Doctor of 
Philosophy in Engineering Mechanics, Virginia Polytechnic Institute, Blacksburg, 
Virginia, June 1968, 



SECTION I. INTRODUCTION 


The problem of determining the natural frequencies of the free vibra- 
tion of a rotating, beam-connected space station is considered by developing a 
mathematical model of the system which represents the general three- 
dimensional motion of the various components of the space station. The con- 
figuration studied is composed of two space modules connected by a flexible 
beam where the system is made to spin in the plane of its orbit in order to 
produce an artificial gravity environment within the space modules. The 
orbital configuration of the space station is shown in Figure 1. 



FIGURE 1. ORBITAL CONFIGURATION OF THE SPACE STATION 

The kinetic energy and potential energy of the space station are 
formulated in terms of a set of generalized coordinates, and a Lagrangian 
function is developed for the system. Hamilton's principle is applied to 
determine the governing equations for the motion of the rotating space station. 
Boundary conditions representing the clamped-clamped attachment of the beam 
to each space module are applied to the ends of the beam. Thus the mixed 
problem of a continuous beam with two large end masses is reduced to the 
problem of a continuous beam with nonhomogeneous boundary conditions. 
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The vibration of the space station is considered to be limited to small angular 
and linear displacements from the motion corresponding to steady rotation as 
a rigid body. Within the limits of this small deflection approximation, the 
motion of the space station in the plane of rotation is shown to be uncoupled 
from the motion of the space station out of the plane of rotation. 

An exact solution is obtained for the beam deflection in the plane of 
rotation. This solution is substituted into the appropriate boundary equations, 
and a characteristic determinant is developed. In addition, the behavior of 
the characteristic determinant as the frequency parameter takes on negative 
values is investigated with five special cases in which the form of the exact 
solution for the beam deflection is modified. 

A procedure to solve for the zeros of the characteristic determinant is 
programmed for digital solution on the IBM 7094. Numerical results of the 
analysis for a given space station configuration are presented in Appendix C. 


SECTION 1 1. REVIEW OF THE LITERATURE 


During the last decade a considerable amount of emphasis has been 
placed upon creating an artificial gravity environment within a rotating space 
station. Suddath [1], Kurzhals and Keckler [2], Krause [3], Polstorff [4], 
and others have studied various single body problems, while such authors as 
Chobotov [5], Fowler [6], Pengelley [7], Tai and Loh [8], and Targoff [9] 
have discussed the problem of rotation of cable-connected space stations. 

The concept of compression-member-connected compartments has been 
examined by Tai, Andrew, Loh, and Kamrath [10] in a paper in which the 
stability and response of 13 rotating space station configurations was investi- 
gated. The configurations studied included single-cable-connected compart- 
ments, multiple-cable-connected compartments, and compartments connected 
by compression members to a central hub. 

A recent paper by Liu [11] presents an analysis of two cable-connected 
space stations rotating about an axis normal to their orbital plane. By using a 
concept of concentrated fictitious masses and a Galerkin approach, a solution 
was obtained for the free vibration of the rotating system. 


f 


i 

i 
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SECTION 1 1 1. DEVELOPMENT OF THE MATHEMATICAL MODEL 


Motion of the Stages 


In order to study the motion of the rotating space station, the motion of 
each of the two stages of the station is described independently. The position 
of the center of mass of the ith stage at any time is given by 

R. = x.I + y.J + z.K (i = i, 2) (1) 

as shown in Figure 2. The motion of each stage is considered to have three 
translation components such that 

R. = x.I + y.J + z.K . (i = 1,2) (2) 

1111 


T Ti * r « Body - Fmd Principal *«•» 

2 , 

•jT 



The rotational motion of each stage is represented by the vector sum of 
five independent angular velocities shown in Figure 3. 


u = 9l + tyk + p .e + a .e + y.e (i = 1,2) (3) 

i i| ic. is. 

ii 
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where 6 represents the rate of change of the attitude angle of the space station 
in the plane of the orbit, ip represents the rate of change of the elevation angle 
of the space station from the plane of the orbit, j8 is the rate of pitch of the 


ith stage about its center of mass, is the rate of yaw of the ith stage about 


its center of mass, and y is the rate of roll of the ith stage about its center 


of mass. Transformations from each of the five sets of axes to inertia space 
axes are given by 
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and 


ITmlj = 


™.= 

cosa. , 

l 

since. . 0 

i 


-sina. , 

cosa. , 0 


l 

l 


0 ’ 

0 ’ *_ 

[T 5 ].= 

cosy. , 

0 , -siny^ 


0 

1 , 0 


siny. , 

0 , cosy. 

t T 2l] = 

(cos ip) 

, (cos0sin$) 


( -sin ip) 

, ( cosQcosip ) 


< 0 ) 

, ( -sine ) 

r~ ( cosip) 

> 

(cose sin ip) 


(i 


(i 


( -sin^cos/3^) , ( cosdcosipcosf } . - sinesin/3^) 

(sin^sin^) , ( -cosecos^sin/3. - sinecos/Sj) 


(sine sin 

. . . . , (sinecos^cos/^ + cos0sin0.) 
....,( -sine cos ^sin/3. + cos0cos/3^) 

(1 = 1 , 2 ) 


= 1 , 2 ) 


= 1 , 2 ) 
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[T 4 32l] 


i 


(cos^coso; - sin^cos/3.sino!.) 

i ii 

( -cos^sino;. - sin^cos^cosa!.) 
(sin#s in/3.) 


. . . , /cos0sin^coso!. + cos0cos^cos/3.sina!. 


\ -sin0 sin/3, sinct 
. . , /-cosflsin^sina!. + cos0cos^cos/3.cosa;. 


) 


( -cos0 sinzpsinoij + 
-sinGsin/^coso:. 


J.cosofA 

■ ■) 


. . . , ( -cos0cos^sinj3 i - sin0cos/3.) 


, / sin0sin^cosa + sin0cos^cos/3.sino!.' 


V + cos0sin/3.sina!. 
' i 1 


■ ■) 


•( 


-sinQsin^sino;. + sin0cos^cos/3.cosa 


.• + cos0 sin/3, cosa. 

i l 


, (-sin0 cos i/® in/3. + cos0 cos/3.) 


■) 


(i = 1,2) 


% 



[ T 5432l]j = 


(cosipcosa cosy - sin^cos/3.sino.cosy. - sin^sta/3.stay.), • • 
' i 1 111 ii 

( - cosies ino; - staiicos/3.cosa.) , • • 

' i 11 

(cosJpcoso!.sinY i - stai/>cos/3.stao.stay. + sta^sta/3.cosy.) , .. 


, /cos0sin^coso.cosy. + cos0cosi/>cos/3.sina.cosy. 
/ -sin0sirv3.sinoi.cosy. + cos0cosipsin/3.stay. 

\ +sin0 cos/3 .stay. 


,( -cos0sta^sino. + cos0cosi/>cos/3.coso^ 
-sin0sta/3.coso^ 


, /cos0 sin ^coso. stay. + cos0 cos ^os/3^sino. stay. 
[ -sin0 sin/3 ^sino. stay. - cos0 cos ips in/3 .cosy. 


V 


-sta0 COS/3 ^COSy^ 


,/sin0sin!/»coso.cosy^ + sin0cosi/'cos/3.sino.cosyj 




+ cos0 sin/3 ^simxcosy. + sin0cos$sin/3.siny. 
- cos0cos/3 i siny i 


,( -sin0sta$sino. + sin0cos^cos/3.coso^ 


V 


+cos0 sin/3 jCoso. 


, /sin0 sin ^coso.. stay. + sin0cos$cos/3.stao..siny. 

I +cos0sin/3.stao.stay. - sin0cos^sin/3.cosy. 

l i i ii 


V 


+cos0cos/3^cosy. 


(i = 1,2) 


The components of the angular velocity about a set of body-fixed principal axes 
of each stage (e , e , e ) are determined from 

1 . Ct G . 

i i i 
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(4) 


co . = 01 + ipk + fi.e + a .e + y.e = co , e + o> e + co _ e„ 
i l 4 l c. 'i s. i. 1. 2. 2. 3. 3 

1 1 xi xi xx 

(i = 1,2) 


Usxng the transformatxon equations to write k, e,e,e,e,e,e in 

4 c. s. 1.’ 2. 3 

— 11111 

terms of the inertia space axes I, J, K, we solve equation (4) to obtain 


a). = 0 (costpcosa.cosy . - sin^'cos^.sino'.cosv. - sin^sin/3 siny ) 

1. l i ill i i 

(sin^.sina.cosy. - cos/3, siny.) + B.cosa.cosy. -a siny 

ill l l ill i 'i 

(1 = 1,2) (5) 


ox = -0(sin^cosi3.cosy. + cosipsino!.) + ^sin/3.cos ot. - fi.sma. + y 

2. l *i l l l i i *1 

l 

(1 = 1 , 2 ) ( 6 ) 


o)„ = 0 ( sin^sin/3 . cos-y - sin^cos/3 sina.siny. + cosipcosa sin-v ) 

3. l 'i i l 'i i 'i 

(cos/3, cosy^ + sin/S.sino;. siny .) + /3.coso;.siny. + c^cosy. 

(1 = 1,2) (7) 


Configuration of the Space Station 


The position of the space station at a given time is shown in Figure 4, 
where the origin of the inertia space axes is fixed at the center of the earth. 
Points 1 and 2 represent the centers of mass of the stages while points P* and 
P 2 represent the ends of the beam connecting the two stations. Vector di is 
the directed distance from the center of mass of Stage 1 to the connection point 
of the beam and vector d 2 is the directed distance from the center of mass of 
Stage 2 to the point where the beam is connected to Stage 2. The unit vector 
e is along the line from the center of mass of Stage 1 to the center of mass of 


Stage 2. The unit vector e is along the line drawn from P* to P 2 . 
we have PlPz 


Thus 


ii 



Ri + dj = (xt - djcos^sino!! - djsin^cosft cosofj) I 

+( yi - dicos0sin^sino'i + dicos0 cos^cos/3 icoso; j\ J 

- disintfsin^iCOSG!! 

+/ z\ - d 1 sin0sin^sinQ'i+ d^ivBcosipcosfiiCosaii 
l+ d 1 cos0sirv3icoso'i 


IT = R^ + d 2 = (x 2 + d 2 cos^sina; 2 + d 2 sin$cos/3 2 cos o; 2 ) I 

p 2 / w 

+/y 2 + d 2 cos0 sink's ino! 2 - d 2 cos0cos#cos/3 2 cosa! 2 \J 
y+ d 2 sin0 sin/3 2 cosq ! 2 / 

+/ z 2 + d 2 sin0sin^sino! 2 - d 2 s in0 cos ^cos/3 2 coso; 2 \K 
\- d 2 cos0 sin/3 2 cosa; 2 J 


d 12 = d 12 e = R 2 - R t = (x 2 - x t ) I + (y 2 - y t ) J + (z 2 - z A ) K 


I = d e = R -R — d 12 + d 2 - dj 
PlP2 PlP2 PlP2 P2 Pi 

= /x 2 - x t + d 2 cos^sino! 2 + d 2 sin^eos/3 2 cosa; 2 \ I 


+ d^os^sino!! + d^in^cos/StCoso!! 


) 


+ /y 2 - y t + d 2 cos0sin^sino! 2 - d 2 cos0cos^os/3 2 cosa! 2 \ J 
+ d 2 sin0 sin/3 2 cosa! 2 + djCosOsin^sina;! I 

d^osecos^cos/SjCoscei + djSutfsin/SiCosai / 



z 2 - z\ + d 2 sin0sinipsinQ! 2 - d 2 sin0cos^cos/3 2 cosc^\ k” 

- d 2 cos0sin/3 2 coso 2 + dtsinesin^sino!! J 

- dtSinflcos^cos/SiCoso!! - djCosesinSicoso;! / 
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FIGURE 4. POSITION OF THE SPACE STATION 

The unit vectors (e , e , e ) describe the beam orientation. 

_ V PiP 2 w 

Vector e lies in the plane of the orbit and defines the direction of the beam 
w 

deflection w( r, t) , while vector e is orthogonal to e and e and 

v & PiP 2 w 

defines the direction of the beam deflection v(r,t) out of the orbit plane. 
These vectors are derived in Appendix A, from which we write 


e 

v 


I 



x 2 - Xi + d 2 cosi/'sino! 2 + d 2 sin^cos^ 2 coso! 2 
+ dicosi/'sinai + d 1 sini/>cosj3icosa'i 

y 2 - yi + d 2 cos0sin^sino; 2 - d 2 cos0cosi/>cos/3 2 cosa' 2 
+ d 2 sin0sinj3 2 cosa2 + d^osAsin^sinai 
- d^osflcos^cos/Jjcosai + djsinflsin^icosa;! 


I 
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x 2 -xi + d 2 cos^sinc^ + d 2 sin^cos/3 2 coso; 2 
+ djcos^sino!! + diSin^cos/Sjcoso!! 

z 2 - + d 2 sin0sin^sino^ - d 2 sin0cos^cos/3 2 coso! 2 

- d 2 cos0sln)3 2 coSQ! 2 + diSine sini^sino!! 

- disinecos^cos^tcosa!! - dtcostfsin/jicosai 


e 

P1P2 


x 

L 


x 2 - x t + d 2 cos^sinoi 2 + d 2 sin^cos/3 2 cosoi 2 
+ dtcos^sino;! + d^in^cos/jtcoso!! 


J 

+ T 


y 2 - yi + d 2 cos0sin^sinofc - d 2 cos0cos^cos/32cosctfc 
+ d 2 sin0sirv32cosa2 + diCosOsini/fcina'i 
- d 1 cos0cos^cos/3 1 cosa; 1 + diSin^siniSicosQ!! 


+ 


K 

L 


z 2 - Zj + d 2 sin 0 sin^sinQ !2 - d 2 sin0cos^cos/3 2 coso! 2 

- d 2 cos0 sin/32 cosa 2 + dtsinesin^sino!! 

- disin0cosi/'cosi3icoso;i - djcos0 sin/3 icoso!i 


e 

w 


L 


z 2 - Zj + d 2 sin0sin^sino! 2 - d 2 sin0cos^cos/32cosa<2 

- d 2 cos0sin/3 2 coso! 2 + dtsintfsin^sincei 

- djsin0cos^cos/3icosa;i - dicos0sin/3icosa!i 



y 2 - yi + d 2 cos0sin^sino! 2 - d 2 cos0cos</'cos/3 2 coso<2 
+ d 2 sin0sin/3 2 coso!2 + dicos0sin^sino!i 
- d 1 cos0cosi/>cos/3 1 cosa! 1 + diSinQ sin/Sjcoso;! 
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Motion of the Beam 


The position of a point on the beam at any time is given by 


R (r,t) = Ri + di + re +w(r, t) e + v(r, t)e 
C P1P2 w V 


( 8 ) 


where r is the position along the beam longitudinal axes measured from P*, 

w(r. t) is the deflection of the beam in the direction of e measured from 

w 

the vector e* , v(r,t) is the deflection of the beam in the direction of tT 
P 1 P 2 _ v 

measured from the vector e . Using the vector identities of the previous 

article, we write Pl ^ 2 


R = / x t + d t ( -cos^sinai - sin^cos^icoscv! ) 
c 


r 

' + L 


x 2 - Xj + d 2 cos^sino !2 + d 2 sin^cos/3 2 coso! 2 
+ djcos^sino;! + d^in^cos^icoso!! 


+ v' 


+ Yi + di( -cosQ sinipsinai + cosd cosipcosfi t cosa>i\ 
-sin0 sin/3 icosc^ 


+ 


_r 

L 


y 2 - yi + d 2 cos0sin^sina! 2 - d 2 cos0eos^cos/3 2 cosa! 2 
+ d 2 sin0sin/3 2 cosa! 2 + diCosSsin^sincei 
- d 1 cos0cos^cos^{coso; 1 + djsin0 sin/3 icosc^ 


w 

L 


z 2 - z t + d 2 sin0sin^sina 2 - d 2 sin0cos^cos/3 2 coso! 2 

- d 2 cos0sin/3 2 cosa; 2 + d 1 sin0sini/'sino'i 

- d 1 sin0cosi/>cos/3 1 cosa'i - d 1 cos0sin/3icosa' 1 
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V 


L 2 


x 2 - x* + d 2 cos^smar 2 + d 2 sin^cos/3 2 cosa 2 
+ dtcos^sina!! + djSin^cos/3 jcosai 


y 2 - yi + d 2 cos 0 suupsina 2 - d 2 cos0 cos ^cos/3 2 coso' 2 
. + d 2 sin0sin/3 2 cosa 2 + dicos0sm$sino!i 

- djCosecos^cos/SiCosofi + diSinOsin/StCosaii 




z i + di/ -sinesin^sino!! + sinflcos^cos/SjCosa!! 
I + cos0 sin/Stcoso!! 



z 2 -Zj + d 2 sin0sin^sina? 2 - d 2 sin0cos^cos/3 2 coso! 2 

- d 2 cos0sin/? 2 cosa! 2 + d 1 sin0sinipsino' 1 

- djsinflcos^cos^cosa!! - djcosOsin^cos o'! 


w 


y 2 - yi + d 2 eos0sin$sinQ! 2 - d 2 cos0eos$cos/3 2 cosa! 2 

+ d 2 sin0sin/3 2 cosa 2 + djcos# sin^smc^ 

- djCosScos^cos^iCoso;! + djSinflsiuSj cosa?! 


v 


L 2 


x 2 - xj + d 2 cos^sino$ + d 2 sin^cos/3 2 coSQ' 2 
+ diCos^sina;! + dtSin^cos/Sicosa;! 


z 2 - zj + d 2 sin0sin^sino: 2 - d 2 sin0cos$cos/3 2 cosa 2 
• - d 2 cos 0 sin/? 2 coso; 2 + d 1 sin0sin^sina'i 

- d^inQcosi/'coSjS! cosai - djcos© sin/3 j cosa! 



K 


Using small motion approximations for ip, fii , /J 2 , a it q; 2 , v, w and neglecting 

third order and higher terms, the position vector R (r,t) becomes 

c 
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X! - + d,ji - d,-J - 

H - -.) 


+ Gi\ 


+ y% + djcose 


+ at 2 d 2 ^ + v 


(i-i 


\ 


r 


JL ^0 1 6 , 

pi ~ T a i ~ M 


1 2 


- d 2 cos0 


i, .i 


n-i) 


\ 


f 


* „2 1 2 , 

- 7^ - ipa 2 




d 12 cos0-|- - d^isine ( 1 - 


2 

r 


H 


v L ) + dz/32Sin0 L 
|- Y" j(d i2 - d t - d 2 )sin0 - (dj/?! + d 2 /3 2 ) cos0 J 


" "l 5 ^'^ (d i2 - d i - d 2 > + (djOf! + d 2 Q! 2 ) J (d 12 - dj - d 2 ) cos0 

)(■ i) 

- d 2 sin0^i- ^ip 2 - j/3 2 - ^o| - + d 12 sin0^ 

M)- 


+1 z t + dtSin^l - ^ip 2 - ~Pi - jat - ipa t 


K 


+ d^ieos# I 


w 

+ — 
L 


j^( d 12 - 


d 2 /3 2 cos0 — 


d t - d 2 )cos0 + (d\P t + d 2 /3 2 )sin0 




v 

-i? ; 


ip(d i2 - d t - d 2 ) + (d^ + d 2 o; 2 ) (d 12 - d t - d 2 ) sin0 


( 


(9) 
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Thus the vector velocity of a point on the beam is 


R (r,t) = - Md 1 + d 12T -d lT -d 2 - - ^ 12 - + d T " 1 ( 1 


+ d 2 a 2 - + v 


+/ yi + d 1 cos0(-i/^^ - PA - - ipai - a^ip) (l 


- dj0sin0 


H 


f - -Pi - -at - M l - — 


d 2 COS0( -# - P 2 P 2 - - *P®2 - G!^) — 

+ d 2 0 sine^l - - jj8i - joi - - d 12 0sin0^ 

+ d 12 cos0^- - d^^'coss/l - -dj/Sisino/ 1 - 


V 

+ d 2 /3 2 0cos0 — + d 2 /3 2 sni0 — 

(d 12 - d 1 - d 2 )0cos0 + d 12 sin0 + (djjgj + d 2 /3 2 )0sin 0 
+ d 2 jS 2 ) eosfl 


w 

L 


w 


’~L [^ dl2 " dl ' d 2) sin0 -JdiPi + d 2 /3 2 ) cos0 


'] 


‘ 7 ? 


[-^(d 12 - d t - d 2 ) + (djCKt + d 2 Q! 2 )j 

• F - ( d 1 2 - d t - d 2 )0sin0 + d 12 cos0j 

"?[' dl2 ^ "^ dl2 ' dl “ d 2 ^ + ( d l a i + d 2 « 2 )l 

• jj d i 2 - d i - d 2 ) cos 6 

" "L 2_ ["^( dl2 ~ dl " d2) + (dlQ!l + d2Q!2) ] (dl2 " dl " d 2 )cos Bj 
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H** 


Zj + djSinQ ( -ipip - a^cny - ipa^ - arf) ^1 - 

+ d^cose^l - - ^p\ - ^ a} - ^ 

- d 2 sin0( -# - P $2 - a 2 a 2 - <P <*2 - a 2 ip)^ 

( 1 9 I? I2 \ r 
1 ‘ 2 ^ " I^ 2 " 2 “ 2 ' ^ ) L 

• r r • / r \ 

+ d 12 0cos0 — + d 12 sin0 — - dj/^sin©! 1 - — J 

/ r \ ° r • r 

(1 - — J + d 2 /3 2 0sin0 — - d 2J 8 2 cos0 — 


+ d^jcos# 


w 

+ L 


-(di 2 - d t - d 2 )0sin0 + d 12 cos0 
+ (d 1 /3 1 + d 2 /3 2 )0cos0 + (d^! + d 2 /3 2 )sin0 


[‘ 


in 0 ^ 


(d 12 - d t - d 2 ) cos0 + ( dj/Sj + d 2 /? 2 )sin0 
- -^C d 12 - d t - d 2 ) + (d^ + d 2 o; 2 ) 

•E 


(d 12 - d t - d 2 )0cos0 + d 12 sin0j 

d 12 ^ - ip{d i2 - d t - d 2 ) + (d 1 o! 1 + d 2 a 2 ) 


v 

L 2 


0 


a 2 )j 


K 


• |^(d 12 - dj - d 2 ) sin0 
- jy -^(di 2 - d t - d 2 ) + (d^j + d 2 a 2 )J (d 12 - d! - d 2 )smdj 


( 10 ) 


Constraint Equation for the Configuration 


The vector from the center of mass of Stage 1 to the center of mass of 
Stage 2 has been previously identified as d 12 . Consider the plane formed by d 12 
and k as shown in Figure 5. 



FIGURE 5. THE PLANE FORMED BY d 12 , k SHOWING THE PROJECTED 
LENGTHS OF THE CONFIGURATION 

The angle between the beam longitudinal axis and the plane of d 12 , k is 

a , where 


sina — a 


dislnott + d?sinc 


dfQf] + d ? ofa 

L 


and cosa 1 - 2 = 1 - ^-(d^ + d 2 Q! 2 ) 2 
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so the projection of the beam onto the d 12 , k plane is L coso; . L represents 
the deflected length of the beam. ^ 

Similarly, the angle formed by the line of length L cosa with the line 
drawn between points 1 and 2 is j3, where 


sinjs B = d i COSQ! i sia ffi + d^cosa^singy 

L^cosc? 


- jof) + d^ 2 (l - 


D 


1 '~2i7 (dlQ;i + d2 “ 2)2 


di^i + d 2^2 
L 


(neglecting higher order terms) 


and cosp 
Thus we have 


* ” 2 ^ 2 “ 1 " 2L 2 ^ dl ^ x + d ^ 2 ) 


di 2 (t) = djcosoijcos^! + L^cosq;cos/3 + d 2 coso; 2 cos/3 2 


“ d > (‘ - & - 


+ L 


1 - 2L 2 ^ d i Q! i + d 2 a 2 ) 


2L" 


(di^i + d 2 /3 2 ) ' 


i jr 7 9v\ 2 / 3 w\ 2 1 / i , 1 2 \ 

• 2 J [Tr) + [it) dr + d *( 1 - 2^ - py < u > 


which is the constraint equation of the configuration. 
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Kinetic Energy of the Space Station 


The total kinetic energy of the space station is composed of the kinetic 
energies of Stages 1 and 2 plus the kinetic energy of the beam. Therefore, 


T Total ' I M A R * + ^ 


+ + ^(IBJjcoL + (IB) 2 o;2, + 0 (IB) 


2 B 


3 ^ 3 , 


1 -U 

+ — f mR 2 dr 

2 J o c 


(12) 


where and R 2 are given in equation (2) 

a; , co , cl) are given in equations (5), (6) and (7) 

1 . & , *5 . 

ill 

R is given in equation ( 10) . 
c 


Potential Energy of the Space Station 


For an orbiting space station the gravity field of the earth is associated 
with the centripetal acceleration of the mass center of the space station in its 
rotation about the center of the earth. We wish to study only the motion of the 
space station superimposed on the translation of its center of mass and neglect 
the small effect of the gravity gradient. Therefore, we consider that the poten- 
tial energy of the space station is the internal bending energy of the connecting 
beam, written as 


U 



+ 



and for the beam used in the analysis I =1=1. 

w v 


( 13 ) 
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SECTION IV. DERIVATION OF THE GOVERNING EQUATIONS 


Development of the Lagrangian Function 


We now have the kinetic energy and potential energy of the space station 
given in equations (12) and (13) . We write the Lagrangian K = T - V to obtain 

K = + 2 - (IA) 1 wjj + — (IA)2C02 t + 

+ ~MgR2 + "^(10)1^12 + "(rB) 2 «2 2 + ^(IBjscolg 



Hamilton's Principle 

By application of Hamilton's Principle in accordance with The Calculus of 
Variations we obtain, after substituting the constraint equation and neglecting 
second order terms, 


6 f Kdt = 
ti 


0 
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= J \- M a ' x , - M r ( x , - Lift 

h 


fixidt 


\ r ± t A -i ^ 

]-mL xj - + — L) - -d^ + jd 2 a 2 + / -mvdrl 


+ 17 - 


M A ‘y 1 - M^(y 1 - Le cose - L,0sin0) 


dyidt 


-mL 


y< - -^-d 1 0 2 cos0 - ^-diesine + 7j-d 2 0 2 cos0 + ^-d 2 0sin0 
2 2 2 ^ 

- 77L0 2 cos( 9 - ~L0sin0 + ^dj^^sInO - ^-d t £[*0cos0 

2 2 z * 

1 •• 1-9 

- dj/S^cos© - — du%sin0 - — d 2 /3 2 0 sin0 

2 * 

1 • • *• 1 _ • • 

+ — d 2 /3 2 0cos0 + d 2 /3 2 0cos0 + — d 2 /3 2 sin0 

2 £ i 


+ -m(w0 2 sin0 - w0cos0 - 2w0cos0 - wsin0) dr 
0 


+ J - M B (Zi - L0 2 sin0 + L0COS0) 

ti 


A 

1-mL 


z\ - ^•d 1 0 2 sin0+ •jkl 1 0cos0 + ^-d 2 0 2 sin0 - ^-d 2 0eos0 
2 2 2 " 

- ^-L0 2 sin0 + ^-L0cos0 - ^-d^^cos© - ^d^dslne 
2 2 2 2 

1 .. i *9 

- d 1 j3 1 0sm0+ — djjSjCOse .+ -rd 2 j320 cos0 

2 2 

^ . * • • 1 • • 

+ — da8 2 0 s in0+ d 2 /3 2 0sin0 - -d 2 /3 2 cos0 
2 2 


Sz<dt 


-LJ ^ tt , . . 

+ J\ -m( -w0 2 cos0 - w’0sin0 - 2w0sin0 + wcos0) dr 
0 
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tg, J ^ 

f fM B L ( z^cose - yjsine + L0) - ( LA) t ( 0 +'(3 t ) 

- ( IB) j ( 0 +'$ 2 ) 


-mL 


1 j?** 1 -9* • 1 ^ 2 * * 1 * • i # * 

J d i0 + J d 20 + — L 0 - -yjdjSinO + — yjdjsing 

1 •• ~ . 1 ... 1 .. 

- "^yiLsutf - -ytd^jcose + jy^foeos© 


1 , 1 •• 1 •• /s “ r 
+ — Zjdjcosfl - -z t d 2 co89 + —ZiLcosfl 


- TjrM^iSin# + jzid 2 /3 2 sin0 - ^^0+ 

+ ~ g" d i d 2^2 “ ^d 2 l/0 - ~d A d 2 i3 j + — d 2 /3 2 


i /N/ .• 4 »J.. 

+ — djL^j - — d 2 L0 2 


L 

+ J~ m 

0 


-YiV/cosQ - Zjwsin# - dj\v 


a ■ ■) 


, r r 

d 2 w — + Lw— 

1 j Xj 


dr 


/ | m b l (Xi - L ^‘) + (IAJ^eyj + 0 V 1 - 0 2 ^ - 6> 2 <*i) 


/+ (IA) 2^(# + 0Qi - Yl) - (IA) 3(0Yi + fly* + (fy + ip) 
f+ (IB) i(0y 2 + 0y 2 - e 2 tp - 0 2 q! 2 ) + (IB) 2 0(# + ba 2 - y 2 ) ' 
- (IB) 3 (0y 2 + 0y 2 + ct 2 + ip) 
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*^ d i + d l Ij + 3 L ^- *1^1 + iT L ^ + 2^ d l Q ’l^ dl + 

- d 2<*2^ d i + - ^y'i d iCOS0(^ + a t ) 

+ ^-yid 2 cos0(^ + o- 2 ) - -jz^sinfl ( ip + o'!) 

+ “Zid 2 sin0 ( ^ + 0*2) + ‘^’d^0 2 (^ 1 + a^) + ^d 2 0 2 (^ + 05?) 


— djd20 2 (2^ + oil + a 2) + g diL0 2 ( ip + o?j) 


- ~d 2 L0 2 (^ + a**) 


-m 


-(dj + r)v + y 1 cos0v + z'jSinOv + d^ 2 ^-^ “ 


d^ 2 ^-v - L0 2 jv 


dr 




- — mLd^ZiCose - yjsiitf) 


- J^igdiLe 2 + l-mLdjLe 2 + |-mLd 2 0 2 - |-mLd 1 d 2 0 2 


^1 + ^2 


1 + 'l'JPi + ■ ( IA )i(^ + ^l) - jmLdi^i 
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- 1 


^-mLd 1 d 2 i8'2 - mlY -^-d 2 + -^-diL^e - ^-mLd^di + L)^ 2 ^ 

6 \ a 6 / o 

-mLdjdge^a - m-^ /£-(r - L)w - d^vjdr | 

(‘ + + t^>] - 1 




\+ “niLd 2 


|/3 2 + -£/3i| - -mLd 2 /3 2 /(yiCOS0 + Zl sin0) 


_ 

(‘ - + 


A 

L 


- m -f J wdr 


^ « * * ■ • * 

+ — mLd 2 (z 1 cos0 - ytsinfl) - (IB) 1 (0+/3 2 ) 

Ci 


M T3 d 2 L0 2 + ^-mLd 2 L0 2 + jmLdidz© 2 - |-mLd 2 0 2 


6 


i 

-mLd 2 /3 2 + — mLd^/S i 


(‘ + + 1>J - 3 1 

/|d,d 2 + |d 2 Lj 6 + mL^|d,d 2 + |d 2 L + |d^S J ft 

L 

- jrmlA&^p! - f (-rw - d t 0 2 w) dr 


+ mL 


<5/3 2 dt 


- / 

H 


M B d > 


[( 1+ ^)“‘ + ^ 

+ jmLd,jjl+^di + ^<%J 

d L i 

+ m-^- J vdr + — mLdjC^ + c^) 


(y\cos0 + z’iSin©) 


dc^dt 
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/v j /v m 4 , A 

- j"M B d 1 L0 2 + ^-rnLdjli 2 - -mLd 1 d 2 0 2 + -mlxi^ 2 ^ 
^ + “l) ^ + "l’ Q!2 | + (IA) + * 0y i " ^ - ® 2q! 1 ) 


+ (IA) 2 (0 2 i/'+0 2 a 1 -0 7l ) - (IA) 3 (0y 1 + 0y 1 + i/' + a 1 ) 


|-mLd 1 2 Q' 1 + 1-mLdtX! - ^mLd^ (d x + |-L^ + | 

(?“■ * ? L ) 


-mId 1 d 2 Q’ 2 


- mLdj 


Uip + a t ) Q‘ 


/[< 

0 L 


- m-£- J J (r - L)v + djflV + r0 2 v| dr 


’0 2 vj 


t 2 | 

+ J 

*1 


M B d2 |( 1 + “L)“ 2+ "L^ 


(yjcos0 +z 1 sin0) 


- -mLd 2 (^ + (%) 


+ — mLdo 
2 i 


( 1 + + -f“i 


+ m 


/vdr 


o 2 , T J T fl2 J. i-mT.rt.H-fl 2 _ — mid 2 # 2 ] 


J ~ m J 

Mgd 2 L0 z + — mLid 2 L0 2 + -mLd 1 d 2 0 - g 


[ 

[( 1 + 1 ?)“° + 


J 


] a 2 + “^<*1 | + (I®) l (®72 + ^Y 2 " ) 

+ ( IB) 2(0 2 ^ + - 072) ~ (IB) 3(072 + 072 + ^* + ^ 2 ) 


1 0 * * 1 * * 

-mLd 2 02 - — mLd 2 x 1 + mLd 2 ^ 


* (?*> + i L ) 


6o; 2 dt 
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+ |-mLd 1 d 2 a 1 + mLd 2 ^d 1 + | l^(jp + a 2 ) 0 2 


L 

/ 

0 


- m f ( rv + d t 0 2 v + rOV) dr 


- / 

ti 


+ / 

ti 


-(IA) 1 (0 2 Yi + 0* + 9a>x) + (IA) 2 ( 6 $ + 9ai + eip+ dotx - Yj) 


+ ( IA) 3( 9 yi + + 9 cti) 


-(IBM 0 2 y 2 + 0* + 006) + (IB) 2 (0* +006 +0* + 006 - Y 2 ) 


+( IB) 3 ( 0 2 y 2 + 9 * + 0 06) 


( 0 I "[ Mb ^ + + 'l 1 Pl ) + i mL ^] 

1. £yjcos0 + ZjSirtfJ - m^cos© - yjSin©) 


■ + M B L + mL (|di + §1) 


- mfdi + r)0 + mdx^il — 


0 2 8 W 2 - m(w - w0 2 ) 


8r 


•• j 

+ md 2 /3 2 — 


+ md 


M* - ■) 

,(i + j^e'% + mdqfiVi - EI-0 




dt 


r = 0 


6y t dt 


6V2dt 


Swdrdt 
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^2 / r 

/ / - 

t, 0 L 


a 2 v 
"B 8r z 


+ -^-Oj + ^oi 2 ) - f 


h 

L 


) - ? mL ipj 


^'y^os# + ZtSinaJ 

t- |m b L + mL ^|d t + 9 2 - mV - 

i- m^dj + r) - md^ - md 2 cii 2 ■£ 

9 2 ("- $+ ^ 


- - Ip + -^-Q!! + ( d i + r) - EI-|^ 




dt 


<5vdrdt 


(15) 


Now we note that x lt y is z* are independent generalized coordinates, 
so we have three equations associated with the coordinates, written as 


x i = 


M fi L + mL 


( d - 1 ? L ) 


NL 


7 mLdi " 


mLd 2 " 

^f“ 2 


— f Vdr 

M, J 


T 0 


(16) 


yi = 


ML . M 

— — 0 2 cos 9 + ~rz — 

m t m t 


L0sin0 


. -21- f ( W 0 2 sin 9 - w0cos0 - 2w0cos0 - wsin0) dr 

M J 


T 0 
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2Mr 


1X1 L I— * n * • * 0 

(L + 2dj) 9 z cos9 + (L + 2d t )0 sin0 + (d 2 /? 2 - d 1 /3 1 ) 0 2 sin0 
- (d 2 j3 2 - d 1 jS 1 ) 0COS0 + (d^ - d 2 /3 2 ) sine 
|_+ 2(d 1 ^ 1 -d 2J 3 2 )0cos0 


(17) 


ML ML 

• • B • o B •• 

Zi = — — — 0 2 sin0 - — — — 0cos0 
M m M_ 

T T 


m 


M 


T 0 


J (-w0 2 cos 9 - w0sin0 - 2w0sin0 + wcosS) dr 


mL 

2M„ 


(L + 2d t ) 0 2 sin0 - (L + 2d t ) 0cos0 - (d 2 /3 2 - d 1 (3 i )9 z cos9 
- (d 2 /3 2 - dj/Sj) 0sin0 - (dj/Sj - d 2 /3 2 ) cos0 
+ 2(d 1 ^ 1 - d 2 /3 2 ) 0sin0 


(18) 


Combining equations (17) and (18) , we may write 


Zjcos# - y^m# = - 


ML 

B mL /T . nJ v 
+ ( L + 2dj) 


M, 


mL 

2M_ 


2M n 


0 - ^(d 2 /3 2 -d^OO 2 


(d 1 /3 1 - d 2 /3 2 ) + j ( w0 2 - w) dr 


T 0 


(19) 
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and 


zjSinfl + yjcos0 


ML 

B mL /T nJ v 
+ -zrr~ (L + 2d t ) 


M„ 


2M„ 


0 2 + 


mL ' m 

+ ~rz (dj^i - d 2 /3 2 )0 + 




M. 


jfc-Wt-wi 

T 

L 

J (w0 + 2w0) dr 
0 


(20) 

We also have 0, ip, y 1; y 2 as independent generalized coordinates with 
associated equations given by 


Jm^L L + ( IA) j + (IB) j + mL^d 2 + |l 2 + djL^ 
[’ 


+ J M^L + -jmL(L + 2dj)J (zjcos0 - yjSin0) 

' L 1 L 

1 /»*••• r 

—mL(d l fi 1 - d 2 /3 2 ) + m J wdr (yjcose + zjSin©) + m J (d t + r) wdr 
- 0 J 0 


(IA) j + mLdj 


(i d ‘ + i L ) 




( IB) * - mLd 2 


(i d > + l L ) 


$2 = 0 
(21) 


M b L(x 1 - L^) + (IA) ^(©yj + 0Yi - d 2 ip - d 2 ^) 

+ (IA) 2 0(0^+0o;i-yi) - ( IA) 3( 0Yi + 0Yi + $ + cq) 

+ (IB) 1 (0Y2 + ’0Y2 - e 2 */' - 0 2 a 2 ) + (IB) 2 0(0^ + da 2 - Y 2 ) 

• • •• •• • • 

- (IB) 3 ( 0 Y 2 + 0Y 2 + <P + (x 2 ) 
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mL 


ip(dl + d t L + ~L 2 ^ - x^dj + ^-1^ + ^d^dj + -|l^ 

- d 2 a 2 ^d t + |l) + 9 2 4> (|d 2 - |d,d 2 + |d t L - |d 2 L^ 

+ d^ 2 ^^— dj + — - d 2 0 2 o! 2 ^— d! + 


+ [ — inLfdj ( ip + oiy ) - d 2 (ip + ct 2 )] - m 


L \ 

f vdrl (yicos0 + Zjsin0) 

« / 


+ J m(d 1 + r) (v + oV) dr = 0 

0 


( 22 ) 


- (IA) 2 Vi + (IA) 2 9 > + (IAh'eoi + [ -(lA) t + (IA ) 2 + (IA ) 3 ] 9 ( i + a t ) 

+ [ - (IA) t + ( IA ) 3 ] b \ J = 0 

(23) 


- (IB) 2 y 2 + (IB) 2 0 >+ (lB) 2 6a 2 + [ - (IB)! + (IB ) 2 + (IB) S ] 0 Op + a 2 ) 

+ [ - (IB)! + ( IB) 3 ] 0 2 y 2 = 0, 

(24) 


Boundary Conditions of the Beam 


The boundary conditions represent a fixed attachment point at each end 
of the beam. That is, the beam deflection and slope at each end are considered 
to be consistent with the motion of the respective stages. The derivation of 
these boundary conditions is shown in Appendix B, from which we write 

Pi = ~Zr I-daW.^L.t) + ( L + d 2 ) w, ! (0, t) ] (25) 

L 
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$2 = ~ [ ( L + dj) w, j ( L, t) - d*w, t 
L 


0*1 - ^ [ d 2 v >i(L,t) - (L + d 2 )v,j 

L 

“2 = ~ [ - (L + d t ) v, t (L, t) + dfv, 
L 


The Governing 


(0,t)] 

(26) 

(0,t)] 

(27) 

1 (0, t)] 

(28) 

Equations 



By substituting equations (16) , (19), (20), (25), (26), (27), (28) 
into equations (21) , (22), (23), (24) and the remaining terms in equation 
(15) and by neglecting second order terms, we obtain a set of six independent 
coordinates (0, ip, y lt y 2 , w > v ) which yield six governing equations and eight 
associated boundary conditions. The equations are shown below as 


Cj0 + C 2 w, i ( 0, t) + C 3 w,j ( L, t) + 


L .. L . 

C 7 f wdr + m J rwdr 
0 0 


0 


(29) 


- C§ 0 yi - Cg0y 2 - - 


Csd 2 C fi (d 1 + L) 


0% (L,t) 


-C s (d 2 + L) 


L 


+ 


L 


0Vi(O,t) 


- C IQ d 2 

L 


+ 


C ll( d .l + D 
L 


Vq(L,t) 


Cm(d 2 + L) 
L 



L 


V,!(-0,t) 


L 

C 7 0y 1 - C 8 0y 2 - C 12 f vdr 
0 
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L . „ . L . L 

m J rvdr - C ? 0 2 J vdr - m 0 2 J rvdr = 0 

0 0 0 


( 30 ) 


(IA) 2 ‘ Tl - (IA) 2 0> - - IA) 2 d 2 e‘ v>1 (L,t) + ( . I A M L + g2 ) ov, 1 (o,t) 

L L 

+ c 1 e'ip+ ■^ a flv, 1 (L,t) - C ^ L + ^ flv.^o.t) + c 5 e 2 yi = o 


(31) 


(IB) 2 V 2 - (IB ) 2 dip- (r - B ^ i.Q V;l( o,t) + (IB)2 ^ l +d l ) g'v.^L.t) 

L L 

+ C 8 d'ip+ ( 0 ,t) - C 8 ^ + d -^-bv ;1 (L,t) + C 6 b 2 y 2 = 0 

L L 


(32) 


L o L 

2 ~ '9 m o* r m r •• 

mw + EIw, im - mfl w - C 24 0w, u + — — J wdr - J wdr 


T 0 


T 0 


+ (C 7 + mr)0 


„ . T , , m 2 L d 1 d 2 , 

C 20 (L + d 2 ) + 2M + mdjr 


w > l ( 0, t) 


/•*»> 

L 


„ , rrrLd? ,, 

-C 20 d 2 - oV / r ( L + d l) + md 2 r 


w ’i (L, t) 


- ~ [C 12 ( L + d 2 ) + Czgdi] fiV^O.t) 
L 


+ ~ t c i2 d 2 + C 2 3(L + di)) 0V[(L,t) = 0 

L 


(33) 
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and 


~ m 2 r 

mv + EIv, mi - C u e\, n - — f Vdr - (C 12 + mr)ip 


T 0 


(C 7 + mr) d 2 ip - 


„ , , v m 2 Ldjd2 , 

C 20 (L + d 2 ) + 2M 1 L +mdjr 


v,i ( 0 ,t) 


1 

+ 

rw 

L 


_ , m 2 Ld 2 , T , , , 

C 20 d 2 + 2 M~ {L + dl) ' md 2 r 


V,! (L, t) - — l (C 7 + mr) eV, 4 ( 0 , t) 
L 


-£• (C 7 + mr) b\, 4 (L,t) = 0 

L 


( 34 ) 


w( 0 ,t) = 0 
w(L,t) = 0 


( 35 ) 

( 36 ) 


EIw, i 4 ( 0 , t) + D 46 + D 2 w, 1 ( 0 , t) + D 3 w, 1 (L,t) + D 40 V 1 ( 0 ,t) 

L L md L 

+ D 5 0 2 w, 4 (L,t) + D 6 e 2 f wdr + D 7 / wdr+^Li / rwdr = 0 


0 


0 


L 0 


( 37 ) 


EIw, n(L,t) + D 8 0 - D 3 w, 4 ( 0, t) + Di 0 w, i(L,t) - D 5 0 2 w , 1 ( 0 , t) 

• L L » md L •• 

+ D 12 0 2 w,i (L,t) + Di 3 0 2 / wdr + Dj 4 / wdr - f rwdr = 0 

0 0 L 0 


( 38 ) 


v( 0 ,t) = 0 

v(L, t) = 0 


( 39 ) 

( 40 ) 
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EIv, n (0,t) - 


<J ±±J ll c 5 0 7l + ^C e 0y 2 - (L + d2) C 7 ey 2 


+ ~^C 8 0y 2 - D 1S ‘ip - D i6 0 2 ip + D 17 v, j ( 0, t) + D 18 v, j ( L, t) 
L 


+ DjgeV,! (0,t) + D 2O 0V, 1 (L, t) + -^-C 7 0 2 J vdr + -^J-m0 2 f rvdr 

L 0 L 0 


L . L 

+ D 21 J vdr + -^-m f rvdr = 0 


0 


L 0 


(41) 


EIv, n (L,t) --^C 5 0 Yl+ (L ^ d±) -C 6 0y 2 - -fec 7 0y 1 + (L ^ dl) C 8 9y 2 
L L L L 

- D 22 $ - D 23 9 2 ip - Di 8 v , j ( 0, t) - D 25 v, j ( L, t) - DjoflV.j (0,t) 

j L . L L 

- D 27 0 2 v, j ( L, t) - -^-C 7 0 2 J vdr - -^-m0 2 J rvdr + D 28 J vdr 

L 0 L 0 0 


ill 


m 


L 0 


J rvdr = 0 


(42) 


where the C's, C's and D's are combinations of physical constants given 
below 


C t = M fi L + jmL( L + 2d t ) 

C 2 = MgLL + ( IA)j + ( IB)j + mL(dj + ~L 2 + djL) 
C 3 = (IA) X + mLdj( |di + -L) 
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C 4 = (IB)! - mLd 2 ( + |-L) 


= (IAh - (IA) 3 

= (IB)! - (IB)g 

= (IA)i - (IA) 2 - (IA), 

= ( IB)! - ( ib )2 - ( IB) 3 


= (IA) S + (IB) 3 + mL(df + diL + -L 2 ) + ML 


(IA) 3 + -mLd^dj + -L) - 


mLdiCi 


11 mLd 2 Cj 

Ci! : (IB) 3 - mLd 2 ( ~ di + ^"L) + 


C i2 = m dj 


C 13 = di M b A+-^) + £mL(2 + -^ 


C i4 = d l( M B + ^ mL ) ^ 


m b z + mL (i di + i L ) di ( 1 + ^ 


Gig = M L + mL(jdi+ |Lj d 


C i7 = (IA)i + -mLdi 


Cig = _mLdid 2 
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C 19 = -g-mlxi^d! + L) 


C,^(i mL -M T ) 


C21 - 


C 22 = 


C 23 = 


M B ( 1+ ^) + 2 md2 ] d2 

( 2 dl + 3 L )] d2 (* + - l) 


ML + mL 
B 


( IB )i + 7rmLd 2 


C24 - 1x11^2^2^1 + 3 L + gd 2 ^ 


C 25 = md 2 


mL d 1 

2M„ + L ML 


LL 


= ilLdj 
6 2M 

T 


c 27 = (IA)3 + ~mLdi ^dj + |-L^ 


^28 - (IA)3 + ^IllLd 2 


C 29 = (IA)i - (IA) 2 + mLdj^dj + 
C 30 = ( IB) 3 - mLd 2 ^|d 1 + 
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C31 = ( IB) 3 + 3 mLidl 


C32 = (IB)! - (IB) 2 - mLdz^di + 

(i d « + l L ) 


O oq = ML + mLI 


Cl - c 2 - 




C, = C 3 ( L + ^ - C i mLd i ( L + 2d 2 ) - 


M 


2M t L 


Cq = 


Ca = 


- C & + C 1 I P*^2 ( L + 2d t ) + ° 4(L + d l ) 

L 2M t L L 


.(IA)! + (IA) 2 - (IB)! + (IB) 2 - mL (| d l - | d l d 2 + ^ d l L - ^l) 


mLCi . , , . 

- is^ (d ‘ - d2) 


~ _ / t a \ mLdiCi 

C5 = -C3 + (IA) 2 + 2M 


/V , mLdoCi 

C 6 = -C4 + (IB) 2 - 2M 


C 7 = mdi 


mCt 

M„ 
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C 8 


m 2 L 2 d 2 

4M t 


Cit 


m 2 L 2 d 1 d ? 

'• = w “ - 4M_ 


Co = Ci 


Cio = 


CiC 


1W3 


M„ 


lvl 2 T 2 a 2 
m L dj 

4M m 


- C 15 - C 19 


Cn = 


c i c i4 m 2 L 2 d 1 d 2 
M t + 4M t 


- C ie - c 18 


/v 

c 


12 


mdi 


Ci mL 

M t L - 2 M t 


(dj + L) 
L 


C 


13 = 


m 2 L 2 d 2 
4M t - 


/N/ 

C 


CiCiM m 2 L 2 d 2 


14 “ 


M„ 


4M„ 


- C 22 + C 24 


~ mL^Cj. _ 

15 2M t 


27 


C 16 - 


mLd 1 C 1 

2M_ 


- C 


29 


C 17 = 


m 2 L 2 d 2 

4M„ 


'28 


C - C 1 C 13 

^18 “ 




- C 


15 “ ^29 


ClQ = 


mLdoCi 


„ mLdpCi „ 

c *» ' + 032 


m 2 L 2 do 


- C 2 2 - C32 


M t L 2 M t 


C24 = C 33 " 


M_ + — mL 


M_ \ B 2 


Di = ~ [C 5 (L + d 2 ) - (IA) 2 (L + d 2 ) - Cgdt + (IB)^] 

L 


D 2 = — [C 8 (L + d 2 ) 2 - 2C 9 d 1 (L + d 2 ) + C 13 df ] 


D3 — — [ -Cg( L + d 2 ) d 2 + Cg( L + d^L + d 2 L + 2d^d 2 ) - Cj3( L + d^) dj ] 

2 


[C 10 (L + d 2 ) 2 - 2C 11 (L + d 2 )d 1 + C^ 2 ] 


[-C io( L + d 2 ) d 2 + Cjj(L 2 + d^L + d 2 L + 2djd 2 ) - Cj,j(L + dj)d^] 
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1 


D e = [C 1 2 (L + d 2 ) - C 2 5 cij] 

L 


D 7 = — [C 2 o(L + d 2 ) + C 2 gdi] 
L 


D 8 = ~ [C 9 d 2 - (IA) 2 d 2 - Cg(L + dj) + (IB) 2 (L + di)] 
L 


D 10 = — [-C 8 d 2 2 + 2 C 9 (L + d 1 )d 2 - C l3 (L + di) 2 ] 
L 2 


Dj 2 = [ — C jo d§ + 2Cn(L + di)d 2 - Cj4(L + di) 2 ] 

L 2 


Dj 3 = — [Ci 2 d 2 - C25(L + dj)] 
L 


Dj4 = — [C 20 d 2 + C 2 g(L + d 1 )] 

L 


Du = ^ [Ci 5 (L + d2) + Cjgd^] 

L 


Dig — [Ci e ( L + d 2 ) + C 2 odi] 
L 


Du = [ Ci 7 ( L + d 2 ) 2 - 2 C 9 ( L + d 2 ) d t + C 2 id 2 ] 

L 2 


D 18 = — [-C 17 (L + d 2 )d 2 + C 9 (L 2 + d 1 L + d 2 L + 2 djd 2 ) - C 21 ( L + d t ) d^ 
L 2 
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rv oj 


Djg - ^ [C 18 (L 4 d 2 ) - 2(C 1A + C 9 ) (L 4 d 2 ) d* 4 C 22 d 2 ] 


D 


20 


/V /N-/ 


"CiofL 4- d 2 ) d 2 4 (C^ + C g) ( L 2 4 djL 4 d 2 L 4 2djd 2 ) 

-C 22 ( L 4 dj) dj 


■^21 “ ^ [C2o(L 4 d 2 ) 4 C 2 gdj ] 

L 


^22 “ ^ [Ci 5 d 2 4 C 1 g(L4d 1 )] 

L 


^23 - ^ [Cied 2 4 C 20 (L4d 1 )] 

L 




D 2 g - ^ [ C 17 d 2 - 2Cg(L4d^)d 2 4 C 21 (L4d 1 ) ] 


D 27 - ^ [C 18 d 2 - 2( Cjj + Cg) ( L + dj) d 2 + C 2 2 (L + d!) 2 ] 


^28 ~~ t C-2()Ci2 + C 2 6(L + dj)] 

L 


SECTION V. SOLUTION OF THE GOVERNING EQUATIONS FOR 
THE MOTION OF THE SPACE STATION IN THE PLANE OF ROTATION 


The Uncoupled Motion of the Space Station 


The governing equations derived in the preceding chapter are seen to 
be uncoupled into two types of motion. Equations (29) and (33) represent the 
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in-plane motion of the station with boundary conditions given by equations (35) 
through (38). Equations (30) , (31), (32), (34) describe the out-of-plane 
motion of the station with associated boundary conditions given by equations 
(39) through (42) . We may now treat the in-plane motion of the space station 
independently from the out-of-plane motion and determine a solution for the 
bending motion of the beam in the plane of rotation. 


The Beam Deflection Equation 


We assume that 6 , the angle of rotation of the space station in the 
plane of the orbit, is given by 

0 (t) = fit + r ( t) 


where fi is a constant and r(t) is a small quantity representing the time 
dependent perturbation of 6 . Then the time derivatives are 

b = 12 + r (t) 

9 = T ( t) 

Substituting the above identities into equations (24), (33) , (35) through (38) 
and neglecting nonlinear terms, we obtain 


Cjt + C 2 w, j ( 0, t) + C 3 w, i ( L, t) + C 7 f wdr + m J rwdr 

0 0 


= 0 

(43) 
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m 2 fl 2 ^ ~ 2 L 


•• 2 9 mi£ r m r * * 

mw + EIw, im - mSrw - C 2 4 ^tv, lt + — — J wdr - J wdr 


,V -t[ c 


+ (Cj + mr)r - ~ |C 2 o(L + d 2 ) + 


T 0 


m 2 Ld 1 d 2 

2M„ 


M 


+ mdjr 


T 0 


w, ! ( 0, t) 


r**/ 

L 


, m 2 Ld 2 

nd, - — zr—t- 


-C 20 d 2 - ~ M ~* (L + d t ) + md 2 r 


w, j ( L, t) 


- ~ [C 12 ( L + d 2 ) + CM flV 1 (0,t) 
L 


/v 


+ — [0^2 + C-23( L + d 4 )] £l^w, i ( L, t) - 0 

L 


w( 0 ,t) = 0 

w(L, t) = 0 

EIw, (0, t) + DjV + D 2 w, x ( 0, t) + D 3 W, j ( L, t) + j ( 0, t) 

L L . m(J L .. 

+ Ds^V, t (L,t) + D 6 J 2 2 J wdr + D 7 J wdr + ^ 1 J rwdr = 0 

0 0 L 0 


(44) 

(45) 

(46) 


(47) 


EIw, n ( L, t) + D 8 V - D 3 w, t ( 0, t) + D i 0 w, t ( L, t) - DgflV x ( 0 , t) 

L L md L 

+ D 12 fi 2 w, ! ( L, t) + D 13 S 2 2 J wdr + D 14 f wdr - f rwdr = 0 

0 0 L 0 


(48) 
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By solving equation (43) for r and substituting into equation (44), we have the 
governing equation of the beam deflection in the plane of rotation, given by 


m 2 ft 2 L 


mw+ EIw, lltl - mS2Sv - C24S2 w, n + — J wdr 


T 0 


m 2 C7 ~ 

+ — 1 ( C 7 + mr) 


M 


T C< 


L L 

/ • • n • • 

wdr - — ( C 7 + mr) J rwdr 

0 Cj 0 


C 20/T , , m 2 Ld 1 d ? md^r C? 

L + d 2 ) + “ + J — + — 2 -(C 7 + mr) 

L 2 M t L L Ci 


w, t ( 0, t) 


-C ?n d? m 2 Ld 7 . _ . , md 9 r ~ 

™ 2 - L ( L + dj) + 2— + —*(07+ mr) 


2M rri L 

T 


w, ! (L, t) 


_[C 12 (L + d 2 ) + C 2 3 d^ ]t2 2 w, i ( 0, t) 

L 


+ ^,fCi 2 d 2 + C 2 3 ( L + d j) ] 2 w, j ( L, t) = 0 

L 


(49) 


With the assumption w(r,t) =rj(t)R(r), the governing equation 
becomes 
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2 2 L 

O 2^,13 / r<- , 0 2 «T? t Wt<\ a. -25 IL y R ( r) (jl - 


11117 R ( r) + EIr)R (r) - mfi 2 77 R(r) - C 24 S 2 2 7 ?R"(r) + 


M, 


T 0 


m 2 C? 


M ~ 

T C< 


L ~ E ^ 

*17 J R(r) dr - mr rj f R(r)dr - tj JrR(r)dr 

0 


L 

/ 

0 


L 

f 

Ci 0 


2 L 

-^-rij f rR(r)dr 


%<L + d 2 ) + + C£l 


2M t L 


'1 


r?R’(0) 


mdi mC; 


ri7R’(0) - 


C 2 (& _ m 2 ] " d 2 . ( l + d t ) + 

L 2M t L Ci 


rjR' (L) 


mdo mC; 


Ci 


rTjR'(L) - — tC 12 (L + d 2 ) + C 2 ^ t ] Q 2 rjR' ( 0) ] 
L 


+ [C 12 d 2 + C 23 (L + dj)] t2 — 0 


(50) 


Dividing by 17 and rearranging terras, we write 


R 

v 


mR(r) - 


m 2 C 
M T C 


f R(r) dr 


mr 


j R(r)dr 


m-^ f rR(r)dr 
C A 0 


2 L 
mr 


f rR(r)dr 

Ci 0 


■ ^ ( L + d 2 ) + 
L 


m 2 Ld 1 d 2 

2 M t L 


C 9 C 


2^7 


Ci 


R' (0) 
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m 2 S2 2 


J R(r)dr + -^[C 12 (L + d 2 ) + C 23 d 1 ] S 2 2 R'( 0 ) 


- ~ [C 12 d 2 + C 23 ( L + dj)] S 2 2 R'(L) - EIR™ (r) + mfl 2 R(r) 

L 

+ C 24 S 2 2 R"(r) ( 51 ) 

We now assume rj ( t) =17 e 1 *^ and take two time derivatives to obtain 
• • 2 ipt 2 

77 (t) = -p r) Q e = — p 17 ( t) . Substituting this value, we remove the time 
dependence from the governing equation and write 


Em IV (r) - C 24 fi 2 R"(r) - m( S 2 2 + p 2 )R(r) 







m 2 fl 2 

ivr 


+ ^ ^ ^ ^ f R(r)dr + mp 2 -^ 1 J rR(r)dr 

T Cb 0 Ci 0 


To simplify the algebra we make two identities, where 
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m 2 J2 2 


m 2 p 2 p 2 C 2 

+ ~^r~ + — — L 

M„ ~ 


L ~ L 

J"R(r)dr + mp 2 ^ 1 - J rR(r)dr 
0 C, 0 

(54) 


Thus equation (52) may be written 


EIR rV (r) - C 24 ft 2 R"(r) - m( 12 2 + p 2 ) R(r) + C,r + C 0 = 0 


wl w2 


An exact solution to this differential equation is 


(55) 


R(r) = Cjsinhajr + c 2 cosha 1 r + c 3 sina 2 r + c 4 cosa 2 r 

C 


C < r 
wl 


w2 


m( n 2 + p 2 ) + m( 12 2 + p 2 ) 


(56) 


where a 


C 24 12 2 +*J C 2 2 4 12 4 + 4EIm( 12 2 + p 2 ) 


2EI 


(57) 


a 2 = 




-C 24 0 2 + J c| 4 S2 4 + 4EIm(12 2 + p z ) 


2 2 X 


2EI 


(58) 


The constants C , and C are now expressed in terms of the arbitrary 
wl w2 

constants c (j = 1 - 4) of the solution by substituting equation (56) into 
equations (53) and (54) . Carrying out this substitution, we obtain 


A l C wl - A 2 C w2 = A 3 c 1 + A 4°2 + A 5C 3 + A 6 c 4 


(59) 
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and 


BiC wl + B 2 C w9 = B 3 Cj + B 4 c 2 + B 5 c 3 + B 6 c 4 


w2 


(60) 


where 


A l - 




L (ft 2 + p 2 ) 


CjLO 2 + p 2 


1 ~ 0 ~ 1 Q 

~ - 2 T _W% T 3 




CjL - — C7I/ L - — mL J L - djCj 
2 3 


- d 2 C 4 


r^> /v rv> 

LC 2 " LC 3 


a 2 - 


mp 2 L 

(a 2 + p 2 ) 


£1 + J± 

m 2 


Aq = 


mp 


I rv r-o 2 1 \ 

— LC 7 (cosha 1 L - i) + — m L^LcoshajL. - — sinhajLj 


a A * a A \ a A 

+ (djCj + LO2) 3 -^ + (d 2 Cj + L C 3) a^cosha^L 


A, = 


mp 


1 ~ ~ 1 ~ / 1. l\ 

— L C 7 sinhaiL + — m L( LsinhaiL - — cosha^ + — ) 
a i 7 1 a i \ 1 a i a t ) 

r*+j r*j 

+ (d^Cj + L C 3) ajSinhajL 


Ac - 


mp 


LC 7 (cosa 2 L - 1 ) + a 


i-m L f- 1 

>2 \a 


sina 2 L - Lcosa 2 L^ 


r*s rs j r^j 


4- (^1^1 LC 2 ) a 2 (^2^1 + L C 3 ) 3-2 < ^^® a 2l J 


A 6 - 


mp 


— — LC 7 sina 2 L + mh( -^-cosa 2 L - ~f~ + Lsina 2 L 


1 

— LU 7 sma 2 L + • , 

a 2 a 2 y a 2 

r^> <-v< fw 

(d 2 C 7 + L C 3 ) a 2 sina 2 L 


a 2 


') 
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m(£2 2 + p 2 ) 


p^L _ mVL^ + P^l (gf + £j) 

L 2M t L C, 


'1 


^ 2 c 1 2 L j. ^ 2 c 23 L (^2 + p 2 ^ 


2M r 


^ p 2 C?L 2 + mp 2 C 7 L 3 
2Cj 3Cj 


(O* + p") 


. 2 2 \ mL y — .9 2 \ p 2 C 7L p 2 C 7 

( O 2 + p 2 ) - — (S2 2 +p 2 ) - ^ 7 - 

T Cjm 


2Ci 


a lP‘ C a> (Ltd;) + a l m2 P^ d I d 2 + aip 2 -^£i 
L 2M m L Cj 


rs*> r^j 


2M t L 


a i» 2 5 i; ., d , a,n 2 C 23 di m 2 (fl 2 + p 2 ) 
L L a ‘ M T 


. P 2 3l . L 

ajCj 


a?Ci 


+ coshajL 


_ a lP *C 2 Q d 2 _ ai m 2 p 2 I^ (L + di) 
L 


2M m L 

T 


a 1 p 2 C 3 C 7 + a i n 2 Ci;d 2 + mLp 2 C 7 


v[i 

/N. 

Cl 

,2: 


L 


ajCi 

* 2/^2 . _ 2 V 


a in 2 C» (L d . m‘<n‘ + p‘) ( p 2 S; 

L a ‘ M T a,C, 
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9 l I I cn o I 



Bfi = 


sina 2 L 


a ?P" C 2 0 d .2. + a 2 m 2 p 2 L d 2(L + dj) 


2 M T L 


a 2 p 2 C 3O 7 a 2 12 2 Cj 2 d 2 naLp 2 C7 


Ci 

i 2 f 


a 2 C t 


_ a 2 n z C 2 8 + P z Cf 

L T &2Cj 


+ cosa 2 L 


mp 2 C 7 mp 2 C 7 


0 ~ 

a|Ci 


a 2 Ci 


Solving equations ( 59 ) and ( 60 ) simultaneously, we write 


C . = A 1 c 1 + A 2 c 2 + A3C3 + A 4 c 4 
w 1 


( 61 ) 


and 


C = B lCl + B 2 c 2 + B3C3 + B 4 c 4 
w ^ 


( 62 ) 


where 


A 2 = 


As — 


M = 


Bh = 


A 3 B 2 

+ 

A 2 B 3 

AjB 2 

- 

A 2 Bi 

a 4 b 2 

+ 

a 2 b 4 

AiB 2 

- 

A 2 Bi 

A 5 B 2 

+ 

a 2 b 5 

A l®2 

- 

A 2 Bj 

A fi B 2 

+ 

A 2 B fi 

AiB 2 

- 

A 2 Bi 

B 3 + 

BA 


B, 
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B, = 


_ B4 + BjAg 

b 2 


B s = 


Bg + 6^3 


B, = 


Bfi + 

B, 


Development of the Characteristic Determinant 


By making the substitutions for 6 = fit + T(t) and w = iq^e ** R(r) 

and substituting r from equation (43) into the boundary condition equations 
given by equations (45) through (48), we have, after algebraic simplification 
and neglecting nonlinear terms, 


R( 0) = 0 

(63) 

R(L) = 0 

(64) 


L L 

EER"( 0) + DjR' ( 0) + D 2 R'(L) + D 3 f R(r)dr - D 4 f rR(r)dr 

0 0 


0 


(65) 


L L 

EIR"( L) + DgR'(O) + D 6 R'(L) + D 7 / R(r) dr + D 8 f rR(r)dr 

0 0 


in/ 

The identities D 4 through D g are given below as 


0 


( 66 ) 


Si 



Cl 


56 




The solution for R(r) given in equations (56), (61), (62) is 
substituted into equations ( 63) through ( 66) to give a set of four simultaneous 
homogeneous equations in the unknown p 2 , written as 

a il c l + a 12 c 2 + a 13°3 + a 14 c 4 = 0 
a 21 c l + a 22 c 2 + a 23 c 3 + a 24 c 4 = 0 
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a 31 c l + a 32 c 2 + a 33 c 3 + a 34 c 4 = 0 
a 41 c l + a 42 c 2 + a 43 c 3 + a 44 c 4 = 0 
from which we obtain an associated characteristic determinant 


D(p 2 ) 


a li 

a 12 

a 13 

a 14 

a 21 

a 22 

a 23 

a 24 

a 3i 

a 32 

a 33 

a 34 

a 41 

a 42 

a 43 

a 44 


The elements a. . of the determinant are defined as 


a il 


B, 


l 14 


m( fi 2 + 

P 2 ) 

/V 

B 2 


m( fi 2 + 

P 2 ) 

b 3 


m( £2 2 + 

P 2 ) 

B4 


m( S2 2 + 

V) 


• u t L Aj 4- B-i 

a 2 i = smhaiL + — — \ *- 

1 1 m(S2 2 + p 2 ) 




1 22 


, T LAo + B 2 
coshajL + /r S 


m(Q i + p e ) 


a 23 = sina 2 L + -~3 


LA> + B-) 


m( £2 + p*) 
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= cosa 2 L + 


LAi + B< 


m(S2 + p) 


Di a i + D 2 a 1 cosha 1 L + (cosha^ - 1) + B + p^ ' 


-H 

a i\ 


, T 1 . , \ A^eD, + 6D ? + 3D 3 L 2 - 2D d L 3 

LcoshajL - — sinhatL 1 + — 1 3 . i * - 2 , a — 

1 a 4 1 / 6m( £2 4 + p^) 


Elaf + D^sinhajL + sinhajL 


„ 1 , , 1 \ B ? L(2D 3 - D 4 L) 

LsinhajL - -cosh^L + - ) + 2 m(n 2 + ^j 


A ; >(6D 1 + 6D; + 3 D 3 L 2 - 2 D 4 L 3 ) 
6m( fi 2 + p 2 ) 


~ ~ Do BqL«(2Dq - DxL) 

D*a 2 + D 2 a 2 cosa 2 L - “aT” (cosa 2 L - 1) + ‘ 2m{Q t +J> $— 


sina 2 L - Lcosa 2 Lj + 

a 2 \ a 2 


A 3 (6Di + 6D? + 3D 3 L 2 - 2D 4 L 3 ) 
6m( S2 2 + p 2 ) 


, ~ D 3 

EIa 2 - D 2 a 2 sina 2 L + — a -sina 2 L 

a 2 


D„ /I „ i , . , \ B 4 L(2D 3 -D 4 L) 

- r- cosa 2 L - - + Lsina 2 L + - * , n 1> ^ A“ 


2m( S2 + p ) 


A 4 ( 6 D< + 6D 2 + 3 D 3 L 2 - 2D 4 L 3 ) 


6 m ( + p z ) 
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) ~ ~ d 7 * 

ElajsinhajL + D 5 aj + DgajCoshajL + — L (cosha 1 L - i) 

a l 




+ -5i (^LcoshajL - — sinhajL) + + ? gL ) 

a. t \ 1 a t 1 / 2m( Sr + p 2 ) 


r^j rs*s 


A t ( 6 D s + 6 D fi ~t~ 3 D7L 2 + 2DgL^) 
6m ( £2 2 + p 2 ) 


2 ~ D 7 

ElaJcoshajL + DgajsinhajL + — “-sinhajL 

a i 


r^J rsj 


+ ^-(LsinhajL - — coshajL + — ) + + *fe L) 

a t \ 1 a t 1 a x ) 2m(£2 2 + p <: ) 


rv_> 


A 2 ( 6 P 5 + 6 Dg + 3D 7 L 2 + 2D S L 3 ) 
+ 6m ( £2 + p ) 


r**> J-J 

-EIa 2 sina 2 L + D 5 a 2 + D 6 a 2 cosa 2 L - — 7 (cosa 2 L - 1) 

a 2 




+ -^a- ( — sina 2 L - Lcosa 2 L ) + + 

a 2 \ a 2 ^ ^ / 2m(Sr + p ) 


/>w* 

A 2 ( 6 Dg -f 6 Dg + 3D 7 L 2 h- 2DgL^) 
, + 6m ( ST + p 2 ) 


-EIa 2 cosa 2 L - D 6 a 2 sina 2 L + — 7 -sina 2 L 

a 2 


5a I _L 


a 2 \ a 2 




cosa 2 L - — + Lsina 2 L ) + 
a 2 


B 4 L( 2D 7 + D a L) 


2m ( Q 2 + p 2 ) 


A 4 ( 6D5 + 6 D 6 + 3D 7 L + 2DgL ) 
6m (fi 2 + p 2 ) 


The characteristic determinant has been programmed for solution 
on the IBM 7094 digital computer. The values of p 2 for which D(p 2 ) = 0 
are the natural frequencies of the in-plane motion of the space station, with 
associated mode shapes. Numerical results of this analysis for a particular 
space station design are presented in Appendix C. 

Special Cases for Values of p 2 < - -Q , 2 

In order to investigate the behavior of the characteristic determinant 
when p 2 takes on negative values of magnitude exceeding S2 2 , we recognize 
five regions where the general solution obtained in the previous article is not 
valid. Discussion of these regions follows. 


r '“'2 4 2 

The Case Where p 2 is Between p 2 = -S2 2 and p 2 = ~ • Within 

this region the governing equation is unchanged from equation (55) , but the 
solution is written more conveniently as 


R(r) = CjSinhair + c 2 cosha 1 r + C 3 Sinha 3 r + C 4 Cosha 3 r 

C 


m( f2 


C r 
wl 
7 


77 


w2 


m( ft 2 + p 2 ) 


(67) 


where a t is given by equation (57) and 


*3 


0 

C24£2 


*\l c| 4 fi 4 + 4EIm(S2 2 
2EI 



( 68 ) 


The previously discussed development procedure (Development of the Character- 
istic Determinant) again leads to a characteristic determinant valid within 
the stated region. Numerical results of the analysis of this case are presented 
in Appendix C. 
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iiiiiiiiii mi mill 


The Case Where p 2 = -£2 2 . For this case the governing equation becomes 
EIR^r) - C 24 fl 2 R"(r) + C wl r + C = 0 

for which the solution is 

C r 3 C r 2 

R(r) = c 1 sinha 4 r + c 2 cosha 4 r + c 3 r + c 4 + - s -^ + — — — 

6C 24 S2 2 2C 24 S2 2 


where a 4 


i?24. Q 

El 



(70) 

(71) 


Following the previously discussed development procedure we obtain a charac- 
teristic determinant valid for the case p 2 = -S2 2 . Numerical results of the 
analysis for this case are presented in Appendix C. 


The Case Where p 2 Takes on Negative Values Larger then 
-C^fl 4 - 4 Elm 12 2 


P 2 = 


4 Elm 


The governing equation for this region is given 


by equation ( 55) , but the solution will now be written as 


R(r) = Cjsinhasrsinagr + c 2 sinha 5 rcosa 6 r + C 3 Cosha 5 rsina 6 r 

C 


c < r 

, wl 

+ c 4 cosha 5 rcosa 6 r + 2 , ^ + 


w2 


m^ + p 2 ) ' m( 12 2 + p 2 ) 


where a, = 


C 24 12 2 

4EI 


J- 


m(12 2 + p 2 ) 


4EI 


(72) 


(73) 


a fl = 


-c 24 fl 

4EI 


J 


-m( 12 2 + p 2 ) 
4EI 


(74) 
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The previously discussed development procedure leads to a characteristic 
determinant valid within the stated region. Numerical results of the analysis 
for this case are presented in Appendix C. 


The Case Where p 2 = 


-2 

-C24& 


- 4EImS2 
4EIm 


2 


is 


The governing equation for this case 


EIR IV (r) - C 24 fi 2 R"(r) + 


4EI 


R(r) 


+ C r + C 
wl w2 


= 0 


(75) 


for which the solution is 

R(r) = c 1 sinha 7 r + 

4EIC .r 
wl 

’ C 2 24 t/ 


c 2 cosha 7 r + c 2 rsinha 7 r + c 4 rcosha 7 r 

4EIC „ 

w2 

C 2 2 4 S2 4 


(76) 


where 


a 7 - 


C^l 

2EI 


(77) 


Following the development procedure discussed previously, we obtain a 

~2 4 2 

characteristic determinant valid for the case p 2 = — 24 —— — 

4EIm 

Numerical results of the analysis for this case are presented in Appendix C. 


The Region Where A 7 B 2 - A 2 B 4 is Identically Zero. This region gives rise to 
a singularity in calculating the value of the characteristic determinant; there- 
fore, the simultaneous solution of equations (59) and (60) must be modified. 
In order to simplify the discussion of this point of singularity, the physical 
constants of the space station design given in Appendix C are used as an 
example. 
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For the given data A 2 = 0 . Therefore, the equation A t B 2 - A 2 Bj = 0 
becomes A 4 B 2 = 0 ; and since B * 0 , we must have Aj = 0 . From the definition 
on page 52 we solve for the p 2 for which A 4 = 0 . Thus 


P 


2 

S 


CjL 




-C 4 Lfi 2 

1 Q ~ ~ 

-mL 3 L - d 1 C 1 


^ 2^1 - LC 2 - LC 3 


which falls within the region discussed in the solution of equations (67) and ( 68 ). 
For this value of p 2 equations (59) and (60) become 


0 = A 3 Ci + A 4 c 2 + A 5 C 3 + A 6 c 4 (78) 

and BgCi + B 4 c 2 + B 5 c 3 + B 6 e 4 + B t C wl 

C w2 = B^ (79) 


The identity of equation (79) is used with the solution for R(r) of equation 
(67) and substituted into the boundary equations given by equations (63) through 
( 66 ) to obtain four simultaneous equations, written as 


wi 

e ll c l + e 12 c 2 + e 13 c 3 + e 14 c 4 + e 15 m (fi 2 + p 2 ) 


wl 


e 21 c l + ®22 c 2 + e 23 c 3 + e 24 c 4 + e 25 m (n 2 + p 2 ) 


wi 


e 3 i c i + e 32 c 2 + e 33 c 3 + e 34 c 4 + e 35 m (g 2 “ p 2 ^ 


J wl 


e 41 c l + e 42 c 2 + e 43 c 3 + e 44 c 4 + e 45 m (ft 2 + p 2 ) 


0 


0 


0 


0 


64 



and a fifth equation is obtained from equation (78) , rewritten as 


wl 


e 5i c l + e 52 c 2 + e 53 c 3 + e 54 c 4 + e 55 " p 2 j~ = 0 


Thus the characteristic determinant for this case becomes 


D <p;) = 


e il 

e 12 

e 13 

e 14 

e 15 

e 21 

e 22 

e 23 

e 24 

e 25 

e 31 

e 32 

e 33 

e 34 

e 35 

e 41 

e 42 

e 43 

e 44 

e 45 

e 5i 

e 52 

e 53 

e 54 

e 55 


where the elements e.. of the determinant are defined as 

ij 

Bj> 

11 m( + p 2 )B 2 

s 


512 m( n 2 + 4 p 2 ) b 2 + 1 


B s 

® 13 ~ m(fi 2 + p 2 )B 2 

S 


e 14 = 


B fi 

m(fi 2 + p 2 )B 2 + 1 
s 


e i5 = 


B 1 

B 2 


e 21 = sinhajL + 


B« 

m(S2 2 + p 2 )B 2 
s 
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e 22 

e 23 

e 24 

e 25 

e 31 


coshajL + 

b 4 

m( S2 2 + p 2 ) B 2 
s 

sinha 4 L + 

b 5 

m(J2 2 + p 2 )B 2 
s 

cosha 4 L + 

Bfi 

m( fi 2 + p 2 ) B 2 
s 




<v ~ , Do , _ _ i . D 3 B 3 L 

D iaj + D 2 a 1 cosha 1 L + (cosha^ - i) + m( QSS a +p 2 )B2 

s 

Da / 1 \ D 4 B 3 L 2 

- -^-^Loosha,L - — „inha,Lj - 2m( 4 n 2 a + p ! )B2 


o D? , , D3B4L 

e 32 = | Elaf + D 2 a 1 sinha 1 L + -^smha^ + m( ^ 4 +p 2 )B2 


- — ^-( Lsinha<L - — cosha t L + — ) 
a l \ a i a ! / 


D 4 B 4 L Z 

■,2 


2 m( S 2 + p ) B 2 
s 


e 3 3 = | £^4 + D 2 a 4 cosha 4 L + (cosha 4 L - 1) + 


- •^ 4 -(Lcosha 4 L - — sinha 4 L ) - — — r^f 
a 4 \ a 4 / 2m(Jr 


D.B.I/ 


2m( O ' 1 + p ) B 2 
s 
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e 34 = 


e 35 = 


e 41 = 


e 42 = 


e 43 = 


D 


EIa 4 + D 2 a 4 sinha 4 L + sinha 4 L + Bf 

S 


1 J-J + , U— 5 

m( 

- 1 Lsinha 4 L - — cosha 4 L + — | - 

4 a 4 a 4 ) 


Pi B 6 L 


L a 4 \ 


2m(£2 + p ) B 2 
s _ 


~ ~ DqL 2 

+ d 2 + -a — 


D,L J 


b 2 


(».-¥) 


9 ~ ~ D? 

ElafsinhajL + D 5 a 4 + DgajeoshajL + — L (cosha 1 L - 1) 


a l 


m(S2 


D 7 BoL Do / 1 \ 

zrb — g". „ + — I LcoshajL - — sinha,L 

& + p‘)B 2 a 4 y 1 a 4 1 j 


DflB 2 L 

■>2 -r 


2m( S2 + p ) B 2 
s 


ElajcoshajL + DcaisinhanL + sinhajL + — 

1 1 b 1 1 a t m(ST + p‘)B 2 


j Lsinha 4 L - — coshajL + — ^ 
a i \ a 4 a ly / 


D 8 B 4 L 4 

,2 , „!T 


2m( fi + p ) B 2 

S 


I ’ /V 


D, 


EIa 4 sinha 4 L + D 5 a 4 + D 6 a 4 cosha 4 L + — L (cosha 4 L - 1) 


D 7 B 5 L 
+ - V2 7 1 ? 


a 4 

1 


/n -, ZlT , + | Lcosha 4 L - — sinha 4 L 

m(f^ + p 4 ) B 2 a 4 ^ a 4 4 / 


s 


d 8 b 5 l 

+ - r” 2 , 2 


2m ( £2 + p 4 )B 2 
s 
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^ D d B L 

i44 = | EIa|cosha 4 L + D 6 a 4 sinha 4 L + -^ L sinha 4 L + — £ ^ 2 - 

4 ' 2 


^ / 

°L Y; 

a 4 \ 


Lsinha 4 L - — cosha 4 L + — )+ 
a 4 


a 4 ) 


DrBrL 2 

2m(b z + p 2 ) B 2 ( 




e 45 = 


Dc + Dr + 


D 7 I/ D fi L d 


+ 


B)L 

B, 


ft * ¥) 


e 51 = A 3 

e 52 = A 4 

e 53 = A 5 

e 54 = A 6 

e 55 = 0 


The characteristic determinant has been programmed for solution on 
the IBM 7094 digital computer. Numerical results of the analysis are presented 
in Appendix C. 


SECTION VI. DISCUSSION 


The results of this analysis of the free vibration of a rotating beam- 
connected space station have shown that the motion of the system can be 
considered to be divided into motion in the plane of rotation which is uncoupled 
from the motion of the system in the plane perpendicular to the plane of 
rotation. For motion in the plane of rotation an exact solution for the beam 
deflection has been obtained, and a set of nonhomogeneous boundary conditions 
representing a fixed-fixed attachment of the beam to the space modules has 
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been used to give a characteristic determinant. Computer programs have 
been developed to solve for the zeros of the characteristic determinant and to 
calculate the natural modes of vibration of the rotating system. 

An exact solution for the beam deflection for motion of the system in 
the plane perpendicular to the plane of rotation can be obtained using similar 
methods, but it was not included in this study because of the length of the 
algebraic forms and computer programs. 

A particular space station made up of two manned space modules 
connected by a flexible beam has been studied to provide an example giving 
numerical results of the analysis. 

The natural frequencies and mode shapes for the six lowest modes have 
been calculated for spin rates varying from 0 to 3.5 rpm. The effect of the spin 
rate has been shown to be an increase in the natural frequencies corresponding 
to an "apparent" increase in beam stiffness due to rotation. 


As discussed in Section V, analyses have been made to investigate 
the behavior of the characteristic determinant when p 2 is less than or equal 
to -fi 2 . The numerical results for these special cases demonstrate the 


existence of two rigid body modes with nonzero frequencies p 2 = and 


~2 4 2 

2 - - 4EImfi , . , 

P = ~ 4EIm , respectively. 


For these frequencies the value of the 


characteristic determinant has been shown to approach zero as p 2 — ■ -Q 2 and 


as p 


~2 4 2 

— C 24 f2 - 4EImfl 


while D(p 2 = -SI*) + 0 and 


D ( 


~2 4 

2 _ -Co - 4EImfi 
P “ 4EIm 


4EIm 

=£ 0. Therefore, for each of these frequencies we have the rigid body mode 
R(r) = 0 since the arbitrary constants c. (j = i - 4) of the solution must 

vanish in order to satisfy the boundary condition equations. 


') 


Frequency sweeps have been conducted for negative values of p 2 , and 
no instabilities for the motion of the rotating system in the plane of the orbit 
have been found. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, May 20, 1968 
964-18-01-00-62 
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APPENDIX A 

DERIVATION OF THE BEAM ORIENTATION VECTORS 


In order to describe the motion of the point on the beam, we identify 

an orthogonal set of unit vectors (e , e , e ) as shown in Figure 4, We 

v P1P2 w 

define the unit vector e from the vector equation 

w 

T x <T = d sinpe 
P1P2 P1P2 w 


= -J 


z 2 - Zj + d 2 sin 0 sin$sina :2 - d 2 sin0cos^cos/? 2 cosa! 2 


- d 2 cos0sin/3 2 cosa 2 + d 1 sin0sin^sino; 1 


- djSinOcos^cos/Sjeos^ - d 1 cos0sin0 1 coscq 


+ K | y 2 - Yi + d 2 cos0sin$sino! 2 - d 2 cos0cos^cos/3 2 coso! 2 
+ d 2 sin0sin)3 2 cosQ; 2 + dtCOSOsin^sinG!! 

- d 1 cos0cos!/>cos/3 1 coSQ! 1 + diSinOsinjSjCoso!! 


where the angle p is shown in Figure A-i. Taking the dot product of the vector 
with itself, we obtain the scalar equation 


(d sinp) 2 
PiP2 


( z 2 - 2 + (y 2 - y t ) 2 + d|sin 2 ^sin 2 o! 2 

+ d 2 cos 2 ;/'cos 2 j32Cos 2 a; 2 + d|sin 2 /3 2 cos 2 Q! 2 
+ dfsin^sin 2 ^ + d 2 cos 2 ^cos 2 /3 i cos 2 a! 1 
+ diSin^cos 2 ^ 


7i 



2(z 2 - Zj) 


d 2 sin0sin^sino! 2 - d 2 sin0cos^cos/3 2 coso! 2 

- d 2 cos0sin/3 2 cosa! 2 + diSinQsin^sina!! 

- djSintfcos^cos/SjCos 

- diCosSsin/S^oso!! 


2(Y2 - Yi) 


d 2 cos0sini/«ina 2 - d 2 cos0cos^cos/3 2 cosa! 2 
+ d 2 sin0sin/3 2 coso! 2 + d^osdsinipsinai 

- djCOSeCOSi/'COS/SiCOSQ'i 

+ d^inesin^^osQ;! 


2d 2 sin^cosipcosj3 2 sino; 2 coSQ; 2 + 2d 1 d 2 sin 2 i/®ina!isina' 2 

2d 1 d 2 sini/'Cos^cos/3 1 cosa! 1 sinQ! 2 

2d 1 d 2 sin^cosi^os/3 2 sina' 1 cosa' 2 

2d 1 d 2 cos 2 ^cos/3 1 cosj3 2 coso; 1 cosQ! 2 

2d 1 d 2 sin/3 1 sin/3 2 coso; 1 coso; 2 

2d 2 sinipcosi/>cos/3 1 sino! 1 coSQ' 1 

K 



FIGURE A-l. THE ANGLE u BETWEEN e AND I 

P1P2 



We now identify the angles \ (k = 1, 2, 3) which the vector d 12 

K 

makes with the inertia space axes shown in Figure A-2a. From Figure A-2bwe 
see that 


7r i 

*1 - 2 + . * 


and 


X 2 — 6 when ip is assumed to be a small angle. 

From the sum of the direction cosines we have, neglecting higher order terms, 



FIGURE A -2. ORIENTATION OF d 12 IN INERTIA SPACE 
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Taking the dot product of d 12 with each of the inertia axes, we write 

d i2 • I = d^cosXj = x 2 - x A - -d i2 *P 

d i2 ’ 7 = d 12 cosX 2 = Y 2 - yi ~ d 12 cos0 

d 12 • K = d 12 cos \3 = z 2 - Zj — d 12 sin0 


Substituting these identities into the equation for (d sin/z) 2 , we have 

P 1 P 2 


(d • sinp) 2 
P 1 P 2 


di2 


+ d 2 i/) 2 sin i! Q ! 2 + d 2 cos^o!2 + d^sin^aq + dfcos^o!! 


2 ,/, 2 , 


+ 2dj 2 


d 2 ^sino ! 2 - d 2 cos^ 2 coSG ; 2 + d^sinaj 

- diCOS/S^OSQ!! 


- 2d 2 ^cosj3 2 sinG! 2 cosa! 2 + 2d 1 d 2 ^ 2 sinG; 1 sina! 2 

- 2d 1 d 2 ^cosj8 1 cosQ' 1 sina 2 - 2d 1 d 2 ^cos/3 2 sina! 1 cosa: 2 

+ 2d 1 d 2 cosQ' 1 cosa: 2 cos(^ 1 - /3 2 ) - 2d 2 ^cos/3 1 sino' 1 coso; 1 
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We assume a u a 2 , jSj, / 3 2 are small angles and neglect higher order terms 
to obtain 


(d sinju) 2 — d?2 ■+■ d 2 + dj - 2dj 2 d2 - ^d^dj + 2d^d2 — (d^2 ~ dj - d2) 
P1P2 


Thus d sin/n = (d 12 - d t - d 2 ) =* L for small deflection theory. 
P1P2 

that a similar procedure can be used to show that d • d 

PlPp PlP? 

for small deflection theory. 


We note 

= d 2 - L 2 
P1P2 


The equations for e and e may now be written as 
w PiP 2 


e 

w 


J_ 

L 


z 2 - z i + d 2 sin0sin^sinQ;2 - d2sin0cos!/>cosi3 2 coso!2 

- d 2 cos0sini32COSQ'2 + djsinesini/'sino'! 

- d^inecosipcos /ijcoso!! - djcosesin/^coso!! 


+ 


K 

L 


y 2 - y t + d 2 cos0sin!/)sinQ!2 - d 2 cos0cosi/)cos/3 2 cosQ!2 
+ d 2 sin0sini32COSQ!2 + djcosflsin^sino!! 

- diCosQcosipcos Picosa-i + djSinesinjS^oso;! 


e = — f x 2 - x i d 2 cosi/ , sina'2 + d 2 sine/'cos / 8 2 coSQi 2 

P1P2 L 

+ diCosi/'sino'! + diSin^cos^cosci!! 


+ 


J_ 

L 


y 2 - yi + d 2 cos0 Sinitsina 2 - d 2 cos0cosi|jcos ^2 COSoi 2 
+ d 2 sin0sin j 82COSQ'2 + d 1 cos0sin^sino; 1 
- djcosecos^cos/sjcoso;! + djSinQsin/SjCosoii 



+ 


K 

L 


z 2 - Zj + d 2 sin0sini/®ina! 2 - d 2 sin£)cos^eos/3 2 cosa! 2 
- d 2 cos0sin£ 2 cosa! 2 + djsinflsinipsina!! 


- djSinflcos^cos/StCoso!! - djcos^sin^jcoso;! 


From the orthogonality of the unit vectors it can be shown by algebraic 
manipulation that the third vector of the triad is given by 


e 

v 



x 2 -x i + d 2 cos!/«ina! 2 + d 2 sin^cos/3 2 cosa 2 
+ djcos^sinoi! + d^ini/'Cos^cosaf! 

y 2 - y 1 + d 2 cos9sinipsino! 2 - d 2 cos#cos^cos/3 2 eosa! 2 
+ d 2 sin9sin/3 2 coso! 2 + djCosSsin^sina:! 

- djcosgcos^cos^icoso;! + d^inflsin^cosa!! 



x 2 - x t + d 2 cos^sino; 2 + d 2 sin^cos/3 2 cosa 2 

+ djcos^sinoi! + djSiniZ-coSiSjCosQ:! 

z 2 - Zj +d 2 sin(?sin^sino! 2 - d 2 sin0cosipcosi3 2 cosQ! 2 

- d 2 cos9sin/3 2 coso! 2 + djsinflsinipsina;! 

- disintfcos^cos/jjcosa!! - djCos^sin/SjCoscei 
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APPENDIX B 

DERIVATION OF THE BOUNDARY CONDITIONS 


The relationship between the slope of the beam at the end points and 
the rotation angles of the respective stages is established for the fixed-fixed 
boundary conditions by the following derivation: 

The orthogonal set of unit vectors e , e , e has been described 

v Pip 2 w 

in Appendix A. From Figure B-l we identify the angles ju , p , M made 

_ _ __ V i r l W 1 

by the vector e 0 with e , e , e respectively, where 
2 , v p^p 2 w 


e^ • e 
2 i v 


cosp = ( -cos^sinofj - sini/'cos/3 1 cosa: 1 ) 

V 1 


x 2 - X! + d 2 Cos^sino !2 + d 2 sin^cos/3 2 cosQ;2 
+ d 1 cos^sina 1 + djsin^cos/^cosa!! 


( y 2 - yi) - costfsin^sina!! + cosQcos^cos/SjCoso;! 

- sinesin/^coso!! 


'+ (z 2 - z t ) - sin0sini/'sino' 1 + cosflsin/^coso!! 
j + singcos^coSiSjCoSQ;! 


d 2 sin 2 ^sina 1 sina!2 + d 2 sin^cos^cos/?2Sina! 1 cosQ'2 


d 1 sin 2 i/'sin 2 Q! 1 + d 1 sin^cosi/'Cos/3 1 sina' 1 cosa; 1 
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+ 


+ 


d 2 sin!/'cos^cos^ 1 coSQ' 1 sinQ'2 - d 1 sin 2 /3 1 cos 2 o'i 
d 2 s in/? jsin/^cos oleosa; 2 - d 1 cos 2 i/)cos 2 /3 1 cos 2 o; 1 
d 2 cos 2 i/'cos P jcos P 2 cosa iCosq; 2 
d 1 sin^cos^cos/3 1 sino; 1 coso; 1 



cosu = 
r i 

+ in - n) 


L 


(x 2 


- £ - x -^ ( -cosi/'sino'! 

- cosOsini/'sino'i + 

- sinSsin/SjCoscvi 


- sin^cos/^cosa^) 


cos0 cos ipeos (3 jcosa j 


+ < Z 2 - z l) 


- sin0sin^sino!i + sinQcos^cosiSjCosa! 
+ cosesin/SjCOSQ!! 


+ 


L 


d 2 sina' 1 sina '2 - d t 
d 2 COSQ' 1 COSQ' 2 COS(/ 3 1 - /3 2 ) 


e 


e 

w 


cosu 

w 


( y2 - yi) 

L 


- sin0 sin ^sulci'! + cosOsin^^osaj 
+ sin0cos^cosj3 1 cosa' 1 


+ 


[ll ~ z i). 
L 


cosOsini/'sino'! - cos0 cos ^cos /3 jcos a t 
+ sinesin/SjCoso;! 


+ 


+ 


A- 


in^ (sinj3 1 coso!isinQ!2 - sin/^sinctficosa^ ) 


C0Si/'C0SQ!iC0SQ! 2 SinO2 - /3 1 ) 
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FIGURE B-i. ORIENTATION OF WITH THE TRIAD <T ~e e~ 

l v p t p 2 w 

Using small deflection theory and neglecting higher order terms, we write 


cosp 


(‘ + if) - “*( 'if ) 


COS U ca 1 

r l 


COSjU — 

w 


*‘( 1 + £ ) + *{%)■ 
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Now the angle ju made by the projection of e on the plane formed by 

_ VPl 2jL 
eT , e is related to v sl (0,t) with the equation 
v PiP 2 


so 


tan [ix ] 
vr^ 


v, i( 0, t) = 


cos /x 

V 


COSj u 

r l 



(‘ + 'l) - “ ! (^) 

- 

«(*) 


and similarly we have for the other end of the beam 
v,i(L,t) = -q '2 - a i ('l") 

Also the angle /x ^ made by the projection of e^ on the plane formed by 

e , e is related to w,i(0,t) with the equation 
P 1 P 2 w 


so 


tan [ ix ] 

wr_ 


w, jlO, t) = 


cosu 

w 


COS/U 

r 




( 1 + ’l) + ^(l) 
l" 

'■(*) 


w, ^0, t) 
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and similarly for the other end of the beam 
w » i( L, t) = 2 


From the equations for w,i(0,t) and w,i(L, t) we solve for 
and /3 2 , written as 


Pi = ” 7 [- d 2 w, t (L,t) + (L + d 2 )w 5l (0,t)] 

L 


fi 2 = t( L + d i)w, t(L,t) - djW.^O.t)] 

L 


and from the equations for v,i(0,t) and v,j(L,t) we write 


a x = ~ [ d 2 v, j(L,t) - (L + d 2 )v, t ( 0, t)] 

L 


a 2 = — l-(L + d 1 )v fl (L,t) + d lV)1 (0,t)] 
L 
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APPENDIX C 

NUMERICAL RESULTS OF THE ANALYSIS 
The Basic Solution 


A program has been written in Fortran IV for the IBM 7094 digital 
computer to solve for the zeros of the characteristic determinant described 
in Section V. This program utilizes an iteration procedure to conduct a 
frequency sweep for incremental values of p 2 . When the value of the 
determinant changes sign, the increment is progressively decreased to 
converge to the value of p 2 for which D(p 2 ) =* 0. This value of p 2 is listed 
as p? (j = 1, 2, . . . N) , the jth natural frequency of the system, and an 


associated eigenvector representing the mode shape of the jth mode is 
calculated. A subroutine is then used to calculate the value of R.(r) for 

incremental values of r from 0 to L, and the results are plotted to give the 
jth mode shape of the beam. 


The data given in Table C-I are used for input data to the computer 
program. These data represent two manned space modules launched by 
Saturn-type launch vehicles and then connected by a flexible tunnel. The 
tunnel is 4 feet in diameter with thin wall construction of a steel wire grid 
sealed with a soft polymer membrane. The grid provides bending stiffness 
and astronaut protective structure while the membrane serves as a micro- 
meteorite shield and as closure for a shirt-sleeve atmosphere. 

TABLE C-I. PHYSICAL CONSTANTS OF THE SPACE STATION 


Component 

' 

Mass 

Principal Moments of Inertia 
(slug - ft 2 ) 

Associated 

Distance 

EI 

(lb - ft 2 ) 

Stage 1 

M a = 1552.8 slugs 
A 

(IA), = 9(i0 5 ) 

(IA) 2 = 1(10 5 ) 

(IA) 3 = 9(10 5 ) 

dj = 46 ft 

— 

Stage 2 

M = 1552. 8 slugs 
B 

(IB)! = 9(10 5 ) 

(IB) 2 = 1( 10 s ) 

(IB) 3 = 9( 10 s ) 

d 2 = 46 ft 

— 

Beam 

m = . 1242 slug/ft 




L = 47 ft 

. 184( 10 9 ) 
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A spin rate SI varying from 0 to 3. 5 rpm has been studied. The 
maximum SI of 3. 5 rpm is such that the gravity level is approximately 
0. 1 g at the tunnel end (approximately 0. 3 g in the crew quarters) . For each 
spin rate the natural frequencies of the six lowest modes are tabulated in 
Table C-n. Mode shapes for these frequencies are shown in Figures C-l through 
C-6. The effect of the spin rate is seen to be a slight increase in the natural 
frequencies corresponding to an "apparent" increase in stiffness of the beam 
due to rotation (see Hurty and Rubenstein [12] ) . 

In addition, the value of the characteristic determinant has been 
shown to approach zero as p 2 approaches -SI 2 . Since the basic solution 
is not valid for values of p 2 ^ -fi 2 due to changes in the governing equation 
or the solution, as discussed in Section V, we consider the special cases in 
the following sections. 


o o n ? o Co^" iQr ^ - 4Elm £1^ 

p^ Between = -*** and 3 — — — 

4Elm 

As shown in Section V, the solution is given by equation (67) . The 
characteristic determinant for this solution has been programmed, and a 
frequency sweep for incremental values of p 2 within this region has been 
conducted. The results indicate that D(p 2 ) — 0 at both end points and that 
there are no values of p 2 within this region for which D(p 2 ) = 0. However, 
a singularity exists within this region, as discussed in Section V, page 63. 
The value of D(p 2 ) becomes large on both sides of this singularity, which 
is evaluated in the last half of this Appendix. 



-fl2 


As shown in Section V, page 62, the governing equation and solution 
are given by equations (69) and (70) when p 2 = -SI 2 . The characteristic 
determinant for this case has been programmed and the value of D(p 2 = -fi 2 ) 
found to be nonzero. However, the numerical results obtained from the basic 
solution and the numerical results of the analysis for the special case of p 
~2 4 2 

between - SI 2 and ~ < ~' 2 4 ^ ~~ EIm ^ — indicate that there is a natural mode at 

4EIm 

p 2 = _q 2 . Therefore, since D(p 2 = -SI 2 ) # 0, we must have R(r) = 0 for all r, 
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TABLE C-H. NATURAL FREQUENCIES OF THE SYSTEM 


Mode 

»; 

(rad/sec) 2 

n 2 = o 

( rad/sec) 2 

n 2 = .01343 
(rad/sec) 2 

« 2 = . 02686 
(rad/ sec) 2 

fl 2 = . 04029 
(rad/sec) 2 

fl 2 = . 05372 
(rad/sec) 2 

Q 2 = . 06715 
(rad/sec) 2 

a 2 = . 08058 
(rad/sec) 2 

fi 2 = . 09401 
(rad/»ec) 2 

ft 2 = . 10744 
(rad/sec) 2 

G 2 = .12087 
(rad/sec) 2 

G 2 - . 13430 
( rad/sec) 2 

1. 

Pi 

.861857(10') 

. 870433(1 0 1 ) 







. 930449( 10 l ) 


gw 

2. 

D 


. 143631 (10 3 ) 

.143819(10 3 ) 

.143990(10 3 ) 

144169(10 3 ) 







3. 

p5 

. 153237( 10®) 

. 1 53302 (10 6 ) 

.153368(10®) 

. 153433(10®) 

.153499(10®) 

. 153564( 10*) 

,153630(10*) 

. 153695( 10®) 

.153760(10®) 

. 153826( 10®) 

, 153891( 10*) 

4. 

P4 1 

. 115565( 10 r ) 

.115590(10 T ) 

. 115614( 10 : ) 

. U5638<10 7 ) 

.115662U0 7 ) 

. 115687( 10*) 

. 115711( 10*) 

.115735(10 7 ) 

. 115759( 10 7 ) 

. 115784( 10') 

. 115808( 10 7 ) 

5. 

P5 

. 44441 1(10 7 ) 

. 444435 ( 10 7 ) 

. 444487(1 0 7 ) 

. 444539( 10 7 ) 

. 444631( 10 7 ) 

. 444644( 10 7 ) 

. 444699( 10 7 ) 

. 444749(1 0 7 ) 

. 444802( 10 7 ) 

. 444853( 10 7 ) 

. 444906(10*) | 

— 

6. 

Pi 

.1 

. 121290( 10 8 ) 

.121300(10®) 

.121309(10®) 

.121318(10®) 

.121327(10®) 

.121336(10®) 

T. 

.121345(10®) 

.121354(10®) 

.121363(10®) 

.121372(10®) . .121381(10®) 






















Station 


FIGURE 


FIGURE C 


C-3. NORMALIZED MODE SHAPE CORRESPONDING TO THIRD 
NATURAL FREQUENCY OF THE ROTATING SYSTEM 



-4. NORMALIZED MODE SHAPE CORRESPONDING TO FOURTH 
NATURAL FREQUENCY OF THE ROTATING SYSTEM 
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FIGURE C-5. NORMALIZED MODE SHAPE CORRESPONDING TO FIFTH 
NATURAL FREQUENCY OF THE ROTATING SYSTEM 



FIGURE C-6. NORMALIZED MODE SHAPE CORRESPONDING TO SIXTH 
NATURAL FREQUENCY OF THE ROTATING SYSTEM 


which is the rigid body mode with associated frequency p 2 = -£2 2 . Thus we find 
that the rigid body frequency is not the usual p 2 = 0 for a nonrotating body, 
but rather a finite number equal to the negative spin rate. 


With Negative Values Larger than 


O / 

n 2 . "^24^ " 4Eln>& 2 
P ITm 


As shown in Section V, page 62, the solution for p 2 within this region 
is given by equation (72) . The characteristic determinant for this region has 
been programmed, and a frequency sweep for incremental values of p 2 has 

been conducted. The results indicate that D(p 2 ) — 0 as p 2 — 

and that there are no values of p 2 within this region for which 


— ~ 4EImQ 
4 Elm 

D(p 2 ) = 0. 



-C 2 |il 4 - 4Elmfl 2 
4Eim 


The governing equation and solution for this point are given by equations 
(75) and (76) . The characteristic determinant for this case has been 


programmed and the value of D 


/ '"■'2 4 2 

2 = ~ 4EImC2 

. P ” 4EIm 


) 


found to be nonzero. 


D 




However, the numerical results in the second and fourth sections of this Appendix 

^2 4 2 

indicate that there is a natural mode at p 2 = — -4EImf2 _. Therefore, since 

) 4EIm ’ 

+ 0, we must have R(r) =0 for all r. Thus we have 

~2 4 2 

another rigid body mode with associated frequency p 2 = — ^ - 4EImt2 _ ag 
the third section of this Appendix. 1111 


4EIm 


The Singularity Where AjBg - AgBj =0 

As shown in Section V, page 63, a singularity exists in the calculation 

of the characteristic determinant. A computer program has been written to 

study this singularity, and the results indicate that D(p 2 ) + 0. Thus, Since 

s 

the numerical results of the analysis in the neighborhood of the singularity do 
not indicate a natural mode at the singularity, we conclude that p 2 is not a 
natural frequency of the system. b 
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